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Abstract 
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^**^ , Darboux coordinates are constructed on rational coadjoint orbits of the positive 

frequency part g+ of loop algebras. These are given by the values of the spectral parame- 
ters at the divisors corresponding to eigenvector line bundles over the associated spectral 
curves, defined within a given matrix representation. A Liouville generating function is 
obtained in completely separated form and shown, through the Liouville-Arnold integra- 
^ ' tion method, to lead to the Abel map linearization of all Hamiltonian flows induced by 

O^ , the spectral invariants. The results are formulated in terms of sheaves to allow for singu- 

jiz ' larities due to a degenerate spectrum. Serre duality is used to define a natural symplectic 

f— ^ , structure on the space of line bundles of suitable degree over a permissible class of spectral 

curves, and this is shown to be equivalent to the Kostant-Kirillov symplectic structure on 
OQ , rational coadjoint orbits, reduced by the group of constant loops. A similar construction 

CT^ ■ involving a framing at infinity is given for the nonreduced orbits. The general construc- 

tion is given for g = gl(r) or sl{r), with reductions to orbits of subalgebras determined as 
invariant fixed point sets under involutive automorphisms. As illustrative examples, the 
Oh, case g = s[(2), together with its real forms, is shown to reproduce the classical integration 

methods for finite dimensional systems defined on quadrics, with the Liouville generating 
function expressed in hyperellipsoidal coordinates, as well as the quasi-periodic solutions 
of the cubically nonlinear Schrodinger equation. For g = 51(3), the method is applied 
to the computation of quasi-periodic solutions of the two component coupled nonlinear 
5J] ' Schrodinger equation. This case requires a further symplectic constraining procedure in 

5^ I order to deal with singularities in the spectral data at 00. 
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2 ADAMS, HARNAD AND HURTUBISE 

Introduction. 

In a series of recent papers [AHP, AHHl, AHH2] isospectral Hamiltonian flows 
on rational coadjoint orbits of loop algebras arising from the Adler, Kostant, Symes 
(AKS) theorem were studied. A systematic way of representing both finite dimensional 
integrable Hamiltonian systems and quasi-periodic solutions of integrable PDE's within 
this framework was developed using moment map embeddings. A representation of 
such fiows in terms of rational matricial functions of the loop parameter determines an 
invariant spectral curve and associated linear fiows of eigenvector line bundles [vMM, 
AvM]. The Dubrovin-Krichever-Novikov technique [K, KN, Du] may be used to solve 
the equations of motion in terms of theta functions [RS, AHHl], but the Hamiltonian 
content of this approach is not evident. 

The link between the algebro-geometric method of integration and the Hamil- 
tonian point of view was made within the context of differential algebras by Gel'fand 
and Dickey [GD, D], the linearization being based upon the Liouville quadrature 
method. Within the loop algebra setting, however, the symplectic content of the 
algebro-geometric integration method has until now only been developed in specific 
examples, without a general theory encompassing all generic cases. Certain aspects, of 
course, are easy to see; the constants of motion, for example, are essentially the coef- 
ficients of the defining equations of the spectral curves, and their level sets are affine 
subvarieties of the associated Jacobians. The algebro-geometric approach thus gives a 
foliation of the phase space by tori, which may be seen in the symplectic framework as 
Lagrangian. However the Hamiltonian content of the actual integration procedure re- 
quires further clarification within the general framework of loop algebras. (For related 
developments in the case of rapidly decreasing boundary conditions, see [FT, BSl, 
BS2], where the link between the scattering transform and Darboux coordinates on 
Poisson Lie groups is developed.) 

In this paper we use the Lie algebraic and algebro-geometric structures associated 
with isospectral fiows of the AKS type on the dual g"*"* of the positive frequency half 
of a loop algebra g to introduce a set of Darboux coordinates on rational coadjoint 
orbits: the "spectral divisor coordinates" associated with line bundles over spectral 
curves embedded in a specific ruled surface. These essentially give a complete separa- 
tion of variables in the Hamilton-Jacobi problem for all AKS fiows on such orbits - or, 
more precisely, a Liouville generating function for the canonical transformation to lin- 
earizing coordinates in completely separated form. The transformation is expressed in 
terms of abelian integrals and hence this procedure derives the Abel map linearization 
within the Hamiltonian setting provided by loop algebras, with the Liouville-Arnold 
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torus identified with the Jacobi variety of the underlying spectral curve. The general 
approach is developed only for the loop algebras 0l(r)"'", sl(r)"'", but the integrated 
AKS fiows for more general loop algebras may be obtained by restricting to invariant 
symplectic submanifolds obtained as fixed point sets under involutive automorphisms. 

Section 1 gives the classical Hamiltonian approach to the introduction of spectral 
divisor coordinates on rational coadjoint orbits and the construction of the completely 
separated Liouville generating function (eq. (1.76)). Certain partial reductions of the 
orbits that arise in most applications are also dealt with. This section is essentially 
classical in spirit, and may in principle be understood without going beyond the tools 
of nineteenth century mathematics (Riemann surfaces, abelian integrals. Lie Poisson 
structures and Hamilton- Jacobi theory). The only additional notion needed as a unify- 
ing factor is that of finite dimensional Poisson subspaces of the dual of a loop algebra, 
consisting of elements that are rational in the loop parameter. Such spaces can, how- 
ever, be given a natural Poisson structure without any further familiarity with loop 
algebras. 

More specifically, the spectral curve S of A/"(A), an r x r matrix valued rational 
function of the complexified loop parameter A, is determined by the characteristic 
equation 

det(Ar(A) - AC/) = 0, (0.1) 

and so sits naturally in C^ or, after compactification, in an ambient ruled surface T. 
The line bundle corresponding to Af{X) is associated to a divisor on the curve iS, given 
by the finite zeroes of a component of the eigenvector, and hence determines a set of 
functions given by the corresponding coordinate pairs (A^, Cm)a'=i,--5 id ~ genus of S) 
on the ambient surface. (The bundle in question is actually of degree g + r — 1, but the 
associated divisor can be normalized with r — 1 points chosen over A = oo.) The key 
point in understanding the link between the symplectic and algebraic geometry is that 
the set of functions {X^Xfi)n=i,...9 turn out in general to form a Darboux coordinate 
system on a reduced version of the coadjoint orbits (Theorem 1.4), and can very simply 
be augmented using invariants related to the spectrum over A = oo to form a Darboux 
system on the full orbits (Theorem 1.5). From this fact, the canonical linearization of 
flows induced by the spectral invariants leads, through the Liouville-Arnold method 
(eqs. (1.82a,b)), to the Abel map (Theorem 1.6) and hence, to ^-function formulae for 
the integrated flow (Corollary 1.7). 

The development of Section 1 is essentially self-contained and sufficiently explicit 
to lead directly to the examples of Section 3, but a deeper understanding of the under- 
lying approach requires the constructions of Section 2. These are aimed at an intrinsic 
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explanation of the underlying symplectic structure from an algebro-geometric viewpoint 
and involve some more modern machinery, such as Serre duality. It also is convenient 
to formulate the results in terms of sheaves, in order to have a setting in which smooth 
and singular curves may be treated on an equal footing. Such a generalization is re- 
quired to allow for the types of spectra occurring in the rational matrix valued functions 
that arise in some of the more interesting applications, such as the coupled nonlinear 
Schrodinger (CNLS) equation. 

The naturally defined class of spectral divisor coordinates, leads to the following 
intrinsic characterization, identifying the (Lie algebraic) Kostant-Kirillov symplectic 
form on rational coadjoint orbits in terms of purely algebro-geometric data. On the 
reduced orbits, the infinitesimal variations of the spectral curves correspond to sections 
of the normal bundle of 5 C T, constrained to vanish at the poles of A/'(A). The space of 
such sections may be identified with V := H^{S, Ks); i.e. to sections of the canonical 
bundle. On the other hand, variations of the line bundles, representing the tangent 
space to the isospectral foliation, are given by the cohomology group W := H^{S, O), 
and Serre duality tells us that W = V*. The tangent space to the reduced orbit is 
thus identified with V (B V*, which leads to a natural symplectic form corresponding 
to this decomposition. The remarkable fact is that this coincides with the Kostant- 
Kirillov form (Theorem 2.7) on the orbits reduced under the action of the subgroup 
of constant loops. A similar construction, involving supplementary data over A = cxd, 
holds on the nonreduced orbits as well (Theorem 2.8). The AKS theorem, together 
with this identification between Lie algebraic and algebro-geometric symplectic forms, 
seems to be at the root of the ubiquitous presence of algebro-geometric constructions 
in the theory of integrable systems. 

Although the main contents of Sections 1 and 2 concern the loop algebras Ql{r)^ 
or 5l{r)^ , reductions to other algebras, obtained as fixed point sets under involutive 
automorphisms, are also placed in the symplectic framework (Theorems 2.9-2.11). This 
allows a determination of the flows induced by spectral invariants through restriction to 
the corresponding invariant symplectic submanifolds, but does not produce an intrinsic 
formulation in terms of separating Darboux coordinates on the coadjoint orbits of 
these subalgebras. The analogous coordinates have yet to be derived in the general 
case, though one obtains, by restriction from the case 0l(r) + , extensions of the results 
of Section 2 applicable to these subalgebras. 

In Section 3, these results are illustrated in a number of examples. As a first 
application, Darboux coordinates are computed on generic coadjoint orbits of the Lie 
algebras 5[(2, C) and sl(3, C), viewed as rational coadjoint orbits of the corresponding 
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loop algebras having only one simple pole. Next, the sl{2) case with n simple poles is 
shown to reproduce the standard linearization results for well-known classical examples 
of finite dimensional systems (cf. [M]). For this case, the "spectral divisor coordinates" 
are essentially just the hyperellipsoidal coordinates, the reductions corresponding to 
fixing spectral data at oo lead to constraints defined by quadrics, and the spectral 
curve is always hyperelliptic. This case also includes the "finite gap" quasi-periodic 
solutions of familiar sytems of PDE's such as the cubically nonlinear Schrodinger (NLS) 
equation. The higher rank "spectral divisor coordinates" are thus really generalizations 
of hyperellipsoidal coordinates. 

Finally, as an illustration of the s[(3) case, involving trigonal curves, the finite 
gap solutions of the coupled 2-component nonlinear Schrodinger (CNLS) equation are 
also obtained by the Liouville- Arnold integration technique. In this case the particular 
structure of the spectrum at infinity leads to further singularities in the curve, and hence 
a decrease of the arithmetic genus relative to the generic case, and an incomplete set of 
Darboux coordinates on the coadjoint orbit. This example is used to indicate how such 
problems may be dealt with by restricting to an invariant symplectic submanifold on 
which the spectral curves share the same generic type of singularities. The linearization 
on the constrained manifold then proceeds in the same way as in the unconstrained 
case. 

Background and Acknowledgements: For a more complete account of the moment map 
construction leading to isospectral flows on rational coadjoint orbits and the algebro- 
geometric method of integration, the reader should consult [AHP, AHHl, AHH2]. 
The present work is the first complete account of the spectral Darboux coordinate 
construction and the Liouville-Arnold integration method on loop algebras, but earlier 
summaries and announcements of the main results communicated at various conferences 
and workshops may be found in [H, AHH3, AHH4, AHH5]. The authors are pleased 
to acknowledge helpful discussions with L. Dickey, B. Dubrovin, H. Flaschka, P. van 
Moerbeke, E. Previato and A. Reyman relating to this material. 
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1. Darboux Coordinates and Linearization of Flow. 

la. Rational Orbits and Spectral Curves 

The Hamiltonian systems to be considered here involve isospectral flows of matri- 
ces determined by equations of Lax type: 

^ = [BW.A-CAj], (1.1) 

where A/'(A), B{X) are rxr matrices depending on a complex parameter A. The matrix 
A/'(A) is taken to be of the form 

"^ N- 

U{\) = \Y + \Y,Y^^ (1-2) 

1=1 

where Y G 0l(r), and {a^ G C}i=i^...^^ are constants. Thus, we are considering rational 
A/'(A) with fixed, simple poles at the finite points {ai} and possibly at oo. Rational 
matrices with higher order poles may be dealt with similarly, but will not be considered 
here for the sake of notational simplicity. 

The particular form (1.2) arises naturally as the translate by Ay of the image 

i=l 

of a moment map from a symplectic vector space parametrizing rank-r perturbations 
of a fixed N x N matrix with eigenvalues {ai}i=i,...n into the dual (g"*")* of a loop 
algebra, represented by r x r matrix functions of the complexified loop parameter A, 
holomorphic in a suitable domain [AHP, AHH2]. This serves to embed a large class 
of integrable systems as Lax pair fiows in (0^)*. The image space for such maps is a 
Poisson subspace of (g^)*, with respect to the Lie Poisson structure, the symplectic 
leaves (coadjoint orbits) consisting of rational functions of A. Since a specific r x r 
matrix representation is involved, we view g as a subalgebra of 0[(r, C) or 5[(r, C), 
obtained generally by reductions under involutive automorphisms (cf. Section 2c). 

The loop algebra elements X G Ql{r) are viewed as smooth maps X : S^ ^-^ flK^) 
from a fixed circle 5"^ in the complex A - plane, containing the points {ca} in its 
interior, and the subalgebra Ql{r)~^ consists of those ^(A) that extend as holomorphic 
functions to the interior of S^. The loop group Gl{r) similarly consists of smooth 
maps g : S^ \-^ G/(r), while the subgroup Gl{r)~^ consists again of those g{X) that 
extend holomorphically inside 5"^. The subspace 0t(r)_ C 0t(r) of loops extending 



DARBOUX COORDINATES 7 

holomorphically outside 5"^ to oo is identified with a dense subspace of the dual space 
Ql{r)~^* through the dual pairing: 

<^,X >:=-!./ tr(MA)X(A))^, (1.4) 

2m J g\ A 

;UG0[(r)_ , Xegl(r) + . 
The matrix B{X) has the form: 

^(A) = (rf$(Ar(A)))+, (1.5) 

where $ G X(g[(r)*) is an element of the ring of A(i*-invariant polynomials on 0[(r)* 
and the subscript + means projection to the subspace 0[(r) + . In general, no notational 
distinction will be made between 0[(r)"'"* and 0l(r)_. The coadjoint action of Gl{r)~^ 
on rational elements A/q of the form (1.3) is given by: 

g:Ql{r)- — >Ql{r)_ 

^:AV— ^^AV ^^ '\ '^^ " . 1.6 

^ A - cii -^ A - cti 

Equation (1.1) is Hamilton's equation on the coadjoint orbit Q^v;, C 0t(r)"'"*, with 
respect to the orbital (Kostant-Kirillov) symplectic form Worb, corresponding to the 
Hamiltonian: 

(^(^) = $(^ + Ay). (1.7) 

The Poisson commutative ring of such functions on Qj^f^ will be denoted jFy . According 
to the "shifted" version [FRS] of the Adler-Kostant-Symes theorem [A, Ko, S], such 
systems generate commuting Lax pair flows. Moreover, they may be shown to be 
completely integrable on "generic" coadjoint orbits [RS, AHP, AHHl] in gt(r)"'"*. 
On such orbits, the (AKS) ring of commuting invariants is generated by the coefficients 
of the characteristic polynomial of A/'(A). 

In analyzing the spectrum, it is convenient to deal with matricial polynomials in 
A, so we define 

C := ^Ar(A) 
A 

= ya(A) + LoA"-i + --- + L,_i, (1.8) 

n 

a(A):=^(A-«^)• (1.9) 

i=l 



where 
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The matrix 



Lo= limAfo{X) = y^N, (1.10) 



i=l 



may be viewed as a moment map generating the conjugation action of Gl{r) on gt(r) 



+*. 



Gl{r) X 0l(r)+* — > gKO^* 

(<7, X(A))^<7^(A)<7-^ (1.11) 

The matrix C satisfies the same Lax equation (1.1) as A/'(A), and the coefficients of its 
characteristic polynomial: 

P(A,z) :=det(£(A)-z/) (1.12) 

generate the same ring of invariants as that of A/'(A). 
Remark: It is also possible to view 

X-^+^[C-Ya{\)]:=C{\) (1.13) 

directly as an element of an orbit in 0l(r)"'"* (polynomial in A~^). Since the ring of 
invariants is the same, the results are equivalent, with a suitable redefinition of the 
Hamiltonians and parametrization of the spectral curve (cf. [AHP]). We retain our 
present conventions, with A/o(A) viewed as the point in Ql{r)^* undergoing Hamiltonian 
fiow, since these are adapted to examining the particular spectral constraints occurring 
at the finite values {A = a^} that appear in specific examples (cf. Section 3). 

The spectral curve Sq C C^ defined by the characteristic equation 

V{\z) = Q (1.14) 

is invariant under the AKS Hamiltonian fiows. Let m be the degree of C{X), {m = n if 
Yj^Oorm = n — lifY = 0) and let {ki} denote the ranks of the matrices {A^i}i=i,...n 
in (1.2) (coadjoint invariants, and hence invariants of any Hamiltonian fiow in QA/b)- 

Lemma 1.1. The spectral polynomial V{\ z) has the form: 

r 

V{\ z) = i-zY + z^-'V,{\) + Y, MX)V,{X)z^-^, (1.15) 

where 

n 

A,(A) := JJ(A - «,)"^^-(0'^-'=») (1.16) 
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and 



deg Vj{X) = 2^ min{j, ki) — j{n — m) =: Kj (1-17) 



1=1 



^ 


1 




z 


A = 




z = 






A ' 




Xm 



Remark: This means that 'P(A, z), and aU its partial derivatives in A or z up to order 
r — ki — 1 vanish at (a^, 0). 

Proof. This follows immediately by expanding det(£(A) — zl) and using the fact that 

C{ai) = Ni n"=i,j/i("i - «j) has rank ki. D 

The structure of V{X, z) implies that on Sq, z ~ 0(A"^) as A ^ oo. This suggests 
assigning z a homogeneity degree m, thereby giving P(A, z) an overall degree rm. We 
may then compactify Sq, regarding it as the afRne part of an r— sheeted branched 
cover of P^, by embedding it in the total space T of 0{m), the ?7zth power of the 
hyperplane section bundle over P^, whose sections are homogeneous functions of degree 
m (cf. [AHHl] and Sec. 2). The pair (A, z) is viewed as the base and fibre coordinates 
over the affine neighborhood Uq := 7r~^(P^ — {oo}). Over Ui := 7r~^(P^ — {0}), we 
have coordinates (A, z) related to (A, z) on Uq fl Ui by: 

(1.18) 

Re-expressing (1.14) as a polynomial equation in (A, 5) extends Sq to Ui, thereby 
defining its compactification S G T. 

Let us assume that S has no multiple components. Let (Ao,-2o) belong to S, and 
suppose that the multiplicity of the eigenvalue zq of C{Xo) is k > 1. It follows from the 
constructions of [AHHl] (cf. also Section 2) that there is a partial desingularisation S 
of S such that, generically, the number of points (with multiplicity) in S over (Aq, zq) 
equals the number of Jordan blocks of C{Xo) with eigenvalue zq, and S is smooth over 
(Ao,-2o)- If, for example, £(Ao) has only one Jordan block of size k with eigenvalue 
Zq, then S = S and, generically, S has a smooth /c-fold branch point over P^. In 
the opposite extreme, if £(Ao) has k independent eigenvectors with eigenvalue zq then 
generically there are k points (with multiplicity) over (Aq, zq) and S has a k-fold node. 

These remarks are of particular importance when Aq = cti, since the Jordan form 
of C{ai) is an invariant of the coadjoint orbit. Thus, if C{ai) is diagonalisable with 
multiple eigenvalues, the generic spectral curve for the orbit will be singular. 

Genericity Conditions 

In what follows, we only consider the singularities that follow from the specific 
structure (1.2), (1.8) assumed for £(Ao). We shall make the simplifying assumption 
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that the A^^ (and hence £(««))) are diagonahzable, with the only multiple eigenvalue 
being z = 0, with multiplicity r — ki. This property is, of course, "generic" for orbits 
with rank(A^i) = ki, but is only assumed in order to simplify the exposition. If other 
Jordan forms are allowed for the A^i's, the only effect is to change the specific form 
(1.22) for the spectral polynomial ^(A, z), (1.27) for the genus formula determining the 
dimension of Q^vb and the explicit expressions (1.83), (1.85) for the abelian differentials. 
All these can easily be modified to hold for other cases. The main results, contained 
in Theorems 1.3-1.6, Corollary 1.7 and the subsequent sections, remain valid mutatis 
mutandis. 

We also assume that one of the following two conditions hold: 

Case (i): Y = Q and Lq has a simple spectrum (to = n — 1). 

Case (a): Y j^ and has a simple spectrum {m = n). 

Again, these conditions are generic and invariant on coadjoint orbits, but in section 3 
it will be indicated how they may be relaxed. 

Finally, we make a further spectral genericity assumption regarding the singu- 
larities of the curve S; namely, that the only singularities occur at the points (0^,0), 
where there is an r — /c^-fold node with r — ki distinct branches intersecting transversally. 
This amounts to requiring that the eigenspaces of C{X) all be 1-dimensional except at 
X = ai, where, by the structure of A/o(A), the eigenvalue z = has an eigenspace 
of dimension r — ki. The desingularization S is then smooth and is isomorphic to S 
away from (a^, 0). This condition is generic in the space of A/q's of the form (1.3) and, 
if satisfied at any point of QA/b? it is also valid in a neighborhood of the isospectral 
manifold through that point. (In particular, it is invariant under the AKS fiows.) 

The coefficients of the polynomials Vj{X) generate the AKS ring on each coadjoint 
orbit Qtv;, C 0l(r)+* and should be viewed as functions on the Poisson submanifold 
consisting of rational elements of the form (1.3) (with rank(A^i) = ki). Note that 

Pi(A)=tri:(A), (1.19) 

and hence its coefficients are Casimir invariants (i.e. constants on all coadjoint orbits). 
The nonzero eigenvalues {zi^^} i=i,...n,K=i,...ki over the points {A = ai}i=i,...n are also 
Casimir invariants, since they are determined as the nonzero roots of the characteristic 
equation: 

n 

det[N, Yl (a^-aj) -zl]=0, (1.20) 
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which is invariant under the coadjoint action (1.6). The A^ := Yl7=i ^« trivial invariants 
{ziK,} determine, in particular, the coefficients of Pi (A), since 

^z,, = tr£(a,)=Pi(a,), i = l,...n. (1.21) 

(For case (i), this is sufficient to determine the degree n — 1 polynomial "Pi (A); for case 
(ii), the degree n coefficient is just tr Y.) This may all be summarized by noting that 
the spectral curves S on the orbit Qf/^^ are constrained to pass through the N + n points 
{(tti, ZiK,), (tti, 0)}, with r — ki branches intersecting at the singular points {(cti, 0)}, the 
values {zi^} being fixed. It should also be noted that for y 7^ the leading (deg Ki) 
terms in the polynomials 'Pj{X) are constants, determined entirely by the symmetric 
invariants of y. For Y = 0, the leading terms are not constants, but they are determined 
as symmetric invariants of Lq, and hence are constant on its level sets. 

A way to express V{X, z) in terms of independent, non-Casimir invariants is to 
choose a reference point Mr G Qa^o 01^ the orbit and parametrize the difference between 
P(A, z) and its value Vr{\, z) at Mr. 

Proposition 1.2. In a neighborhood of the point Mr G Qmo' ^^e characteristic poly- 
nomial has the form: 

r 

ViX,z)^VRiX,z) + aiX)Y,cijiX)PjiX)z'-^ (1.22) 

i=2 



where 



a,(A) = Y[{X - «^)™^-(o,J-fc.-i)^ (1.23) 

p,{X)=:Y,PjaX'' (1-24) 



a=0 

and {Pj(A)}j=i,...r are polynomials of degree: 

. . ( d,~j tf Y = 
6, = deg Pi(X) = { ^ 1.25a) 

' ^^'^ ' \dj tf Y^O ^ ^ 

n 

dj = y^ min{j — 1, hi). (1.25b) 

For y = 0, the leading coefficients Pjs^ are constant translates of the elemen- 
tary symmetric invariants of Lq, while for Y j^ 0, the leading coefficients Pjs are 
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all constants; namely, the elementary symmetric invariants of Y (translated by the 
corresponding leading terms in Vr{\,z)). The number of spectral parameters {Pja}, 
fa = 0, . . . 5j •+ n — TO — 1, j = 2, . . .r) defining the polynomials pj (A) on generic orbits 
is thus: 

r 

d = ^ {dj - {n - m){j - 1)) 
= g + r-l, (1.26) 



where 



1 1 " 

g=-{r- l)(mr - 2) - - ^^(r - k,){r -h- 1). (1.27) 

In a neighborhood of any generic point on Qj^„, these spectral invariants are all inde- 
pendent. 

Proof. The structure of V{X) follows Lemma 1.1, plus the fact that V{X, z) — Vr{\ z) 
vanishes at each A = a^, while z vanishes at least linearly in A — ct^ along each branch 
through (ai,0). From formula (1.6) and the above genericity conditions regarding the 
residues A^^, the dimension of the coadjoint orbit Qj^^ is easily computed to be 2d. From 
the proof of complete integrability of the AKS flows on such orbits given in [AHHl] , it 
follows that the isospectral foliation is Lagrangian, and hence the d spectral parameters 
{Pja} are independent. The expression of Pjs. in terms of the elementary symmetric 
invariants of Lq or Y follows directly from the fact that the leading term in C{X) in 
eq. (1.12) is either LqX^-^ or YY". D 

It follows from the adjunction formula applied to the curve S obtained by blowing 
up T once at each point (ct^, 0) (cf. [AHHl, GH]) that 'g in eq. (1.27) is also equal to 
the (arithmetic) genus of S. If we reduce such an orbit under the G/(r, C) action (1.9) 
for case (i), or the action of the stabilizer Gy C Gl{r, C) of y for case (ii), the dimension 
of the reduced space is precisely 2'g, and the projected spectral invariants again define 
completely integrable Hamiltonian systems [AHHl]. These facts suggest exploiting the 
orbital symplectic structure further so as to explicitly integrate the isospectral flows 
via Hamiltonian methods. This will be the content of the following subsections. 

lb. Divisor Coordinates on Reduced Orbits 

Define 

1C{\ z) := t{\) - zl , (1.28) 
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and let /C(A, z) denote its classical adjoint (matrix of cofactors). Let Vq G C be an 
eigenvector of Lq in case (i), or of Y in case (ii). From the results of [AHHl], it follows 
that the set of polynomial equations: 

JC{X,z)Vo = (1.29) 

have, away from (ct^, 0), precisely 5? generically distinct finite solutions {(A^, -2^)}^=i,...g 
that may be viewed as functions on the coadjoint orbit QaTq- (Changing to the coordi- 
nates (A, 5), there are also r — 1 further solutions with A = 0, i.e., A = 00. If Vq is not 
chosen as an eigenvector of Lq or Y, the remaining r — 1 solutions will generically also 
be at finite values of (A, z).) 

The significance of these functions in terms of the algebraic geometry of the spec- 
tral curves S may be summarized as follows (cf. [AHHl] and Section 2a below for 
the detailed construction). To each matricial polynomial C{X) is associated a degree 
g + r — 1 line bundle E ^ S over the partly desingularized spectral curve S. Away from 
the degenerate eigenvalues this coincides with the dual of the bundle of eigenvectors of 
£^(A) over S. At a smooth point (A, z) of 5, the fibre of E is the cokernel of the map 
/C(A,2): 

— , C'^ ^SIM, C'^ ^ ^ — , 0. (1.30) 

More generally, this exact sequence defines the direct image of E over S (cf. Sec. 2a). 
Vectors Vq in C^ then give sections of E by projection. These sections vanish precisely 
at the points where Vq is in the image of IC{X,z). Since /C(A, 2)/C(A, 2) = V{X,z)I, 
this is equivalent to (1.29), at least over the open set of points in S corresponding to 
nondegenerate eigenvalues, for which the corank of /C(A, z) is one. From [AHHl], the 
degree of i? is ^r^ + r — 1, so sections of E have g + r — 1 zeroes. The choice of Vq as an 
eigenvector of the leading term in jC{X) implies that r — 1 of these are over A = 00, and 
the coordinates of the remaining g points are the finite solutions {(A^, -2^)}^=i,..,g. 

In evaluating Poisson brackets, it is preferable to introduce another normalization, 
corresponding to the eigenvalues of —y^ rather than £(A), by defining: 

(■■=4vs (1-31) 

and 

-M(A,C):=^-a, (1.32) 

with classical adjoint Ai{XX)- Then 

JC{X,z) = [a{X)Y-'MiXX) (1.33) 
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and eq. (1.29) is equivalent to: 

MiXX)Vo = 0. (1.34) 

The g solutions {{X^, z^)}^=i^,,,g are thus related to the solutions {{X^XiJ.)}n=i,...g of 
(1.34) by: 

Viewing {{X^Xn)}^j.=i,...g ^.s functions on Qa/o? ^^ '^^Y evaluate their Poisson 
brackets with respect to the orbital (Kostant-Kirillov) symplectic structure Uorb- 

Theorem 1.3. The Poisson brackets of the functions {X^Xn)tJ.=i,...g o'^c-' 

{A^,A,} = 0, {C^,G} = 0, {A^,C.} = 5^.. (1.36) 

Proof. Choose a basis in which the leading term in jC,{X) (i.e. Lq for case (i) and Y 
for case (ii)) is diagonal, and let Vq = (1,0... 0)"^. Let A4ij{XX) denote the ij^'^ 
component of A4{X, C). The points (Aj^, Ci/) are then determined by the conditions 

Mki{X.X.)=0 (1.37) 

for all k. Generically, these points are cut out by only two of these equations, say 

Mil = M21 = 0. (1.37a) 

That is, generically the matrix 

dX dC / 

is invertible. By implicit differentiation, the Poisson brackets of the functions (Aj^, Ci,) 
are then: 

{Ajy, A^l {XiyXfJ,} 

.1 /{A^ii(A^,C^), A^ii(A^,Cm)} {-Mii(A^,0, A^2i(A^,Cm)} \ /zr ^T-l 



^ ''' \{M2i{X,X.). MiiiX^X^)} {-M2i(A.,C.), -M2i(A^,Cm)}; (1.39) 

To determine the brackets in the matrix on the right hand side of equation (1.39) 
we first recall that if / and g are functions on the orbit Q^/b 5 their Poisson bracket at 
a point n G Qa/q C 0t(r)_ is given by 

{F,G}=</x,[^,0]>, (1.40) 
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where -^ is the differential of / at /U, considered as an element of g[(r)"'", and the pairing 
< , > is defined by eq. (1.4). The ij^^ coefficient of M.{\,C,) evaluated at the point 
(Ao, Co), viewed as a function of /U G 0l(r)_, may be written 

-M,,(Ao, Co) = - < li. Y^ > - Co/, (1.41) 

A — Ao 



where Cji is the matrix with a 1 in the jilh. place and zeroes elsewhere. It follows that 

(1.42) 



5Xij(Ao,Co) _ e^i 



b\i A — Ao 

and hence, dropping the subscripts, 

1 



{{M^l{\, C) - M^l{a, ri))6jk - {Mkj{X: C) - Mkj{(T: V))S^l] ■ 



-^"^ (1.43) 

Since A^(A, C) is the classical adjoint of 7W(A, C) we have 

M{\ 0M{\ C) = det(A^(A, 0)1. (1.44) 

Differentiating with respect to a parameter t yields 

dM{\X) _ ■M(A,C)tr((|.M(A,C)).M(A,C)) -■M(A,C)(iA^(A,C)).M(A.C) 
dt det(7W(A,C)) ■ 

away from points (A, C) where det(7W(A,C)) = 0; i.e., points on the spectral curve. 
Thus, away from the spectral curve, 

^M^J{\0 ^ ■M,p(A, C).M,,(A, C) - M,p{\. C)-M,,(A, Q , . 

dMp,{\X) detAl(A,C) ■ ^ ■ ^ 

The derivation property of the bracket 

{■M.,(A,C),.MM(a,r^)} = ^ ft^^^a ^! l^''^^''^\ {M,,{\X).Mrs{a.v)} (1-47) 

then gives 

1 



{Mn{X,0,Mki{a,7j)} = { 



[{Mia,v)M{XX))kiMn{a,v) 



A — (J [detM{a,r]) 

- {M{a,v)MiXX)hiMki{cr,v)] 



iM{XX)Mia,v))kiMn{XX)] 



(1.48) 
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By eq. (1.37), M.ki{^uXu) vanishes for all k,u. Taking the limits (A, C) -^ (A^jCm)? 
(a,?]) -^ i^iyXu) along any path transversal to the curve S, the right hand side of 
equation (1.48) has limit zero for u ^ n (the simple zero in detA^ is cancelled by a 
double zero in the numerator), implying 

= {M2i{X,, G), M2i{X^, C^)} = (1.49) 

when /U 7^ z/. Hence {A,^, A^}, {C;/, Cm) ^^^ {-^j^' Cm) ^^^ vanish when z/ 7^ (U. 

To compute the bracket for u = n we first note that the brackets on the diagonal 
of the matrix on the right hand side of equation (1.39) are zero in this case. Thus, to 
show that {Xi,Xu} = 1 it suffices to show that 

{Mii{X,Xu).M2i{X,Xu)} = det{F,). (1.50) 

To compute the left hand side of (1.50) we first take the limit (A, () -^ (a, 77) in (1.43) 
using the derivation property (1.47) of the bracket to show 

prs ^ ^ 

On the other hand 

Equation (1.50) now follows by substituting eq. (1.46) into eqs. (1.51) and (1.52) and 
using the fact that A4{X^, (,y) has rank 1. D 

The implication of Theorem 1.3 is that the functions {(A^, Cfi)}fi=i,...g nearly pro- 
vide a Darboux coordinate system on the coadjoint orbit Qj^f^ . However, the dimensions 
are not quite right. For case (i), we have 

dim Q^„ =25^+(r + 2)(r-l) (1.53a) 

for generic orbits, while for case (ii) 

dimQ^,^2{g + r-l). (1.53b) 

(Note that in these formulae, it is the value of ^r^ that is different, according to eq. (1.27), 
not the dimension of Q_\f^^ which, of course, is the same.) 
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On the other hand, for case (i), the Marsden-Weinstein reduced coadjoint orbit 
Qred, obtained by fixing the value of the Gl{r) moment map Lq and quotienting by 
its stabihzer Glo C Gl{r), is of dimension 2g. Similarly, for case (ii) we may reduce 
by the stabilizer Gy C Gl{r) of Y, since the shifted AKS Hamiltonians of the form 
(1.7) are invariant under this subgroup and the restriction of Lq to the corresponding 
subalgebra Qy is conserved under the fiows. The reduced orbit under this action, also 
denoted Qred, is again of dimension 2g. (Note again that the value of g for the latter 
case is, by eq. (1.27), ^r{r — 1) greater than for the former.) Thus, if the coordinates 
(A^,Cm) could be shown to be projectable to the reduced spaces, and if the reduced 
Poisson brackets remain the same as in eq. (1.36), we would have Darboux coordinates 

on Qred- 

For case (ii) this may be seen immediately. Since Vq was assumed to be an eigen- 
vector of y, with no degeneracy allowed, the defining equation (1.34) is invariant under 
the stabilizer Gy C Gl{r) (an r — 1 dimensional abelian group under our hypotheses). 
Thus {Xn,Cn)n=i,...g ctre all invariant under the Hamiltonian Gy- action, and the Pois- 
son brackets of their projection to Qred are the same as on Q^^ . 

For case (i), we cannot quite apply Hamiltonian symmetry reduction under Gl{r), 
since the functions (A^, (^) are only invariant under the stabilizer subgroup Glq- How- 
ever, we may still compute the Poisson brackets on the reduced space by the procedure 
used for constrained Hamiltonian systems. Let us first choose the reduction condition 
given by the level set: 

Lo = diag{/J, (1.54) 

where the eigenvalues {Ii} are, by our genericity assumption, distinct. The diago- 
nal terms in eq. (1.54) are the first class constraints, which generate the Hamiltonian 
Gl^ - action, and the terms with i > 1 may be chosen as the independent generators. 
Applying the standard procedure of modifying the Hamiltonian by adding a linear com- 
bination of the remaining, second class constraints, we see that the following modified 
functions generate fiows that are tangential to the constrained submanifold: 

A, = A^- J2 %37^M,., (1.55a) 

c=c- i: %%^M,.. (1.55b) 

Evaluating their Poisson brackets, we find, again: 

{A^,AJ = 0, {C^,CJ = 0, {A^,C}=5^. (1.56) 
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since, by implicit differentiation of the defining equations (1.37), the second factor 
in (1.55a,b) involves terms of the form {7Wfci(A, C), (-^o)ij} which, applying the chain 
rule and eq. (1.46), vanish unless i = 1. The cross terms in the Poisson brackets 
(1.56) therefore all contain terms proportional to {7Wfci(A, C), (-^0)11}, ^ 7^ 1, which 
vanish at (A, C) = (A^,Cm)- Since the functions (A^,^^) coincide with (A^,^^) on the 
constrained manifold and generate tangential flow, it follows that the projections of 
(A^,Cm) to Qred (the quotient of the constrained manifold by Glq) satisfy the same 
Poisson bracket relations as (1.56). Finally, for other values of Lq than (1.54), we just 
repeat the same argument with respect to a diagonalizing basis of eigenvectors. 

Combining these results we obtain, for both cases (i) and (ii): 

Theorem 1.4. The projections of {\^, C,n)iJ.=i,...g ^o the reduced orbit Qred? in both case 
(i) (Y = Q) and case (ii) (Y ^ Q, with distinct eigenvalues), are Darboux coordinates; 
that is, the reduced symplectic form is: 

9 

t^red = ^ dX^ A c/Cm- (1-57) 

Remark^: The proof of Theorem 1.3 did not depend on the fact that there are g 
finite points in the spectral divisor. If the vector Vq is not chosen as an eigenvector of 
y, the number of such finite points, and corresponding coordinate pairs (A^, C^), may 
be between g and g + r — 1. The number of points over A = 00 equals the number 
of eigenvalues z of the asymptotic form of C{X) (i.e., Y for case (ii) and Lq for case 
(i)), for which Vq is in the image of Y — zl for case (ii) (resp. Lq — zl for case (i)). 
This is zero for generically chosen (non-diagonal) Y (or Lq) or, equivalently, if Y is 
taken as a diagonal matrix, and Vq chosen as a vector with no vanishing components 
(e.g. Vq = (1,1,... 1)^). In this case, the number of finite spectral divisor coordinate 
pairs (A^, C^) will actually be g+r—1, sufficient to provide a Darboux coordinate system 
for the full orbit in case (ii) and an r{r — 1) co- dimensional symplectic submanifold in 
case (i) (cf. Section Ic). However, for the examples involving integrable systems that 
will be of interest to us (cf. Section 3), it is not this type of spectral Darboux system 
that is needed for directly determining solutions, but those derived in the following 
subsection. The problem lies with the invertibility of the Abel map (cf. Section Id), 
which requires a degree 'g divisor. The remaining r — 1 points of the spectral divisor 
are related to the singular differentials having pole singularities over A = 00. 



^ Thanks are due to B. Dubrovin for raising the point discussed in this remark. 
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There remains then the question of the nonreduced orbits Qmo ■ Can the functions 
(A^, Cm)m=1v5 somehow be completed to provide a Darboux coordinate system on Qa^q? 
The answer is: yes, for case (ii), and partially for case (i). The construction is given in 
the following subsection. 

Ic. Darboux Coordinates on Unreduced Orbits 

In case (i) we shall obtain Darboux coordinates, not on the complete coadjoint 
orbit Qa/'o7 but on a constrained submanifold Q^^ C Qj\f^ consisting of elements for 
which the off-diagonal elements of Lq vanish: 

(Lo).j=0 if i^j. (1.58) 

By our earlier genericity assumptions, the diagonal elements (-Lo)ii are hence distinct, 
and it is easily verified that Q^ C Qtv;, is a symplectic submanifold of dimension 

dimQ0^^=2(5^+r-l). (1.59) 

(Note that m = n — 1 in the genus formula (1.27) and we are dealing with case (1.53a), 
not (1.53b).) For case (ii), we choose a basis in which Y is diagonal: 

Y = diag{yj. (1.60) 

Thus in both cases, the leading term of t{\) is diagonal. As in the proof of 
Theorem 1.3, we also choose the eigenvector Vq in (1.29) to be Vq = (I7 0, 0, ..., 0)^. In 
both cases, let 

P, :=(Lo)^^, z = l,...r. (1.61) 

These generate the action of the group D of diagonal matrices, which equals Glq and 
Gy, respectively, for cases (i) and (ii). The generator Pi is not independent of the 
others, since the sum: 

r 

$^P^ = trLo (1.62) 

j=i 

is a Casimir. These generators Poisson commute amongst themselves and also with the 
D-invariant functions {\^.X^l)^l=l,...g^ since equation (1.37), which determines them, is 
D-invariant. In case (i), let 

g, := ln(Li),i + ^ Y, ln(P^-P;•), (1-63) 
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while for case (ii), let 

Q, :=ln(Lo)a. (1.64) 

With these definitions, we have: 

Theorem 1.5. The coordinate functions {X^, Cm' 1i^ Pi)fi=i,...g;i=2,...r form a Darboux 
system on Q^ in case (i), and Qmo, in case (ii); that is, the only nonvanishing Poisson 
brackets between them are given by: 

{A^, 0} = ^tJiu: {q^: Pj} = Sij. (1.65) 

Equivalently, 

9 r 

uJorb = ^ dX^ AdC^ + '^ dqi A dPi, (1.66) 

where the equality refers to the full orbit Q_a/o ^'^ case (ii), and the restriction of Uorh 
to Q^ in case (i). 

Proof. The proof proceeds in two steps. First, as in Theorem 1.3, the Poisson brackets 
are computed on the full coadjoint orbits. In case (i), we then reduce this to the 
constrained submanifold, which is symplectic. From the Poisson brackets (1.43) used 
in the proof of Theorem 1.3 follows: 

{{Lo)ij, {Ls)ki} = {Ls)kj^ii — {Ls)iidjk, (1.67a) 

{M,j{X, 0, {Lo)ki} = {Y- M{X, C))u6jk - {Y - M{X, C))kj6u (1.67b) 

{M,j{X, 0, {Li)ki} = (A - 5^ a^)[{Y - M{X, QhSjk - {Y - M{X, 0)kjSu] 

m 

+ [{Lo)iiSkj - {Lo)kj]Sii. (1.67c) 

This implies, in addition to the relations (1.36), the brackets: 

{A^,?^} = {C^,?^} = (1.68a) 

{X^,P,} = {C^,,P^} = (1.68b) 

{P,,Pj} = (1.68c) 

{q,,Pj} = 5,, (1.68d) 

r (P,-P.)-^ for case (i) 

[ U tor case (nj, 

where (1.68a) holds only on the constrained manifold Q^ for case (i). As in the 
proof of Theorem 1.3, we must use the fact that Ain is zero at (A^,C/j)- In case (ii) 
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this completes the proof. For case (i), the constraints must be taken into account. 
As in the proof of case (i) of Theorem 1.4, we shift the functions (A^, Ci^, Q'i, -Pj) by 
terms proportional to the second class constraints (-^o)ij = 0, i ^ j to get functions 
(A^, Qv, Qij Pj) which agree with (A^, Ci/, Qi, Pj) on Q^ and which generate flows in Q^vb 
that are tangential to Q^ . Since 

{(Lo)ij, {Lo)ki} = {Lo)kj5ii - {Lo)ii5jk, (1.69) 

it suffices, for a general function / on Qj^^ , to take 

/ = /- E %^(^o),.. (1.70) 

As in the proof of Theorem 1.4, the Poisson brackets (1.36) remain unchanged on the 
constrained manifold. From eq. (1.69), it follows (as in Theorem 1.4), that the Pj's 
already generate tangential flows and hence the Poisson brackets (1.68b-d) remain un- 
changed. Eq. (1.68a) also is unchanged since, by the same arguments as in the proof 
of Theorem 1.4, the additional cross terms obtained after constraining are all propor- 
tional to terms of the form {gj, (Lo)ji}, which vanish on the constrained manifold. 
Using eq. (1.67a), we see that the remaining Poisson bracket (1.68e) gets shifted to 
zero. n 

Remarks: 

i) The submanifold Q^ C Qa/^ is, in fact, the relevant phase space for many in- 
teresting examples of integrable systems, such as the flnite gap solutions of the 
cubically nonlinear Schrodinger equation (cf. [AHP, P] and Sections 3c, 3d). 

ii) If, in formulae (1.2), (1.3), we choose n = 1, ai = and 1" 7^ 0, then Qjvb is 
really a coadjoint orbit in gi{r)* or si{r)* and Theorem 1.5, together with the 
Ad* invariants (Casimirs), provides Darboux coordinate systems for these flnite 
dimensional Lie algebras (cf. Sec. 3a). 

Id. Liouville - Arnold Integration and the Abel Map 

We now turn to the integration of the Hamiltonian systems (1.1) generated either 
by elements of the Poisson commutative ring J-'y of functions of the form (1.7), with $ in 
the ring /(^[(r) ) of Ad* - invariants on Ql{r) , or its extension jFy(P) by the generators 
{Pi}i=2,...r- Thus, our Hamiltonians are all expressible as functions of the invariants 
{Pia, Pi}- The notational conventions of the preceding sections allow us to treat cases 
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(i) and (ii) simultaneously, although it should be remembered that the spectral curves 
and ring of invariants J-'y{P) depend on the choice of Y, and the relevant symplectic 
manifold is Q^^ for case (i) and the entire orbit Q^/'o ^or case (ii). The reduced spaces, 
though both denoted Qred, are also different, their dimensions 2g being given by the 
genus formula (1.27) with m = n — 1 for case (i) and m = n for case (ii). For case 
(i), Qred signifies the generic Gl{r) - reduction of Q^v;, or, equivalently, the reduction 
of Q^^j by the abelian r — 1 - dimensional group action generated by {Pi}i=2,...r- For 
case (ii) , Qred is the reduction of the full orbit Qj^^ by the latter action. 

The g + r ~ 1 independent spectral invariants for case (i) may be chosen to be 
{Piaij Pj)i,j=2,...r;a,=o...Si-i, siuce the Coefficients Pjs^ occurring in Proposition 1.2 may 
be expressed as translates of the elementary symmetric invariants of Lq = diagjPi} 

p,s, = i-iy-^ J2 ^^1 • • • ^^. + ^^•' (1-^1) 

where the constants {rrij} depend on the reference polynomial Vr{X,z). For case (ii), 
the leading coefficients {Pjs } are constants (translates of the elementary symmetric 
invariants of Y) and the next to leading coefficients are translates of linear combinations 
of the Pi's: 

r 

P,,s,-i = {-iy-'Y.P^ E Y,,...Y,^_,+n,, (1.72) 

where again, the constants Uj depend on Vr{X,z) and the constants Yi. Thus, the 
g + r — 1 independent invariants may be chosen to be {Pim, Pj}i,j=2,...r,ai=i,...Si-2- The 
Hamiltonians may be viewed in the two cases as functions of the independent invariants: 

h^h{P,a,,Pj) ij = 2,...r, a, -l,...5,-l-e. (1.73) 

with e = for case (i) and e = 1 for case (ii). We can use the Darboux coordinates of 
Theorem 1.5 to expresi 
derivative of a 1-form: 



Theorem 1.5 to express the symplectic forms cUorb or cUorblgo as (minus) the exterior 



e:=J2C^dX^ + J2P^dq^. (1-74) 

Restricting to the invariant Lagrangian manifolds C obtained by fixing the level sets 
of {Pia^Pj}, there exists (within a suitable neighbourhood of such £'s) a Liouville 
generating function 5'(A^, g^, P^a, Pi) such that: 

e\c = dS. (1.75) 
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Integrating from an arbitrary initial point thus gives 

Six,, q,, P,„, P.) = J2 r '^^'^nf'^ dX + J2 Q^P^^ (1-76) 

where the A integrals are evaluated within a chosen polygonization of the spectral curve 
S and the function 

z = z{\P,a,Pj) (1.77) 

is determined implicitly along £ by the spectral equation: 

P(A,^(A,P,„,P,)) = 0. (1.78) 

Applying the standard canonical transformation procedure, the coordinates (QiajQj) 
canonically conjugate to the invariants {Pia, Pj) are then 

dS J^ f^^ 1 dz ,^ , „ , 

fj, — _L /^ 

Evaluating the integrands by implicit differentiation of eq. (1.78) with respect to the 
invariants {Pia^Pi}, and using eqs (1.71), (1.72), we have 



— -ayX) — — r — , z = 2, ...r, a = 1, . . .o^ — 1 — e 



dP^a ' ' V,{\Z) 



(1.80a) 



B. _ ,,,)^ %%(A)(-^)---A--- _ ,^,,..., ,,,„,) 



dPi ^ ^^ V,{\z 



where 



Rij :— < 



(1.81) 



and e = for case (i) 

and e = 1 for case (ii) 
The flow is then given in implicit form by the linear equations: 

J^ f^^^R^^ai(X)(-zY-^X^^-'^ dh ,^ „^, , 
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where {Cja, Ci}i=2,...,r; a=i,...,Si-i-e dtJCB integration constants and a fixed base point Aq 
has been used in the integration. 

Remark: On any given level set of the P^'s the Hamiltonians h{Pia, Pj) project to the 
reduced space Q^ed and eq. (1.82a) alone gives the corresponding linearization of the 
reduced fiow. 

We note that the linearizing map defined by eqs. (1.82a,b) involves 'g + r — 1 
abelian integrals on S. If we knew that the 'g differentials 

a,(A)z^-^A" ^, 
^la ■■= ^ ,. — ^—dX (1.83) 

appearing as integrands in (1.82a) were all holomorphic (i.e. abelian differentials of the 
first kind) and independent then, up to a normalizing change of basis, (1.82a) would 
just be the statement that the Abel map 

A:S~3S — >J{S) (1.84) 

taking the unordered set of g points {p^} G S with coordinates (A^, z^) to its image in 
the Jacobi variety J{S) linearizes the flow - the familiar type of result usually obtained 
from algebro-geometric methods of integration [AvM, KN, Du, AHHl]. In fact, this 
is exactly the case. Moreover the remaining r — 1 integrands 

appearing in (1.82b) are abelian differentials of the third kind with simple poles at the 
r points {ooi} over A = oo with local coordinates 

f (A = 0, 5 = Pi) for case (i) 
oo, ^<^ ^- ' ^' (1.86) 

I (A = 0, 2 = Yi) for case (ii). 

Theorem 1.6. The g differentials {uJia}i=i,...g in eq. (1.83) form a basis for the space 
H^{S,K^) of abelian differentials of the first kind (where K^ denotes the canonical 
bundle). The linear flow equation (1.82a) may therefore be expressed as: 

A(r') = B + Ut, (1.87) 

where B, U G C^ are obtained by applying the inverse of the g xg normalizing matrix 
M, with elements 

M^,(,„) := / u,a. (1.88) 
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to the vectors C,H G C^ with components Cia and —-^p^, respectively (the pair (ia) 
viewed as a single coordinate label in C^ ). 

The r — 1 differentials {uJi}i=2,...r in ^Q- (1-85) are abelian differentials of the third 
kind with simple poles at oo^ and ooi, and residues +1 and —1, respectively. After a 
suitable translation by elements of H^{S, K^) to obtain the standard normalization with 
respect to a canonical homology basis {a^,b^ G Hi{S,Z)}^=i^^^g, these provide a basis 
for the r — 1 dimensional space of normalized differentials with simple poles over A == oo. 

Remark: Combining these results with the remark following eq. (1.82a,b), we see that 
for Hamiltonian flows on the reduced orbit (or equivalently for Hamiltonians that are 
independent of the Pi's), the linearization map only involves abelian differentials of 
the first kind. For flows on the unreduced orbit it is necessary to introduce the differ- 
entials of the third kind in order to determine the time dependence of the additional 
coordinates {qi} {viz. Corollary 1.7). 

Proof. Every holomorphic 1— form on S can be obtained by evaluating the Poincare 
residue of a meromorphic 2— form on T with pole divisor at S. Over the affine coordi- 
nate neighborhood Uq such a residue has the form 

-=|^. (1.89) 

Vz{X,z) 

where, for holomorphicity at A = oo, the total weighted degree of the polynomial 
f{XjZ) must not exceed m{r — 1) — 2, and for holomorphicity at the points {(ai,0)}, 
the function f{X,z) must vanish to sufficiently high order so as to cancel the zeroes 
of Vz{Xj z). Since at these points Vz{Xj z) vanishes like (A — aiY~^\ while z vanishes 
along each intersecting branch like A — ctj, /(A, z) must be a sum of terms of the form 



JJ(A-aj)"^=''^(°'*-'=^-i)A" i = 2. 



where, in order to have total degree at most rM(r — 1) — 2, 

0<a<5j-l-e. (1.90) 

But these are precisely the 1— forms iVia of eq. (1.83), which therefore span the entire 
5'— dimensional space of holomorphic 1— forms H^{S,K^). Since, by Proposition 1.2, 
there are exactly g such w^'s, they are necessarily linearly independent. 

Turning to the remaining r — 1 differentials {ooi} of eq.(1.85), these have the same 
structure near the points (ai,0) as the uJiaS, and hence are holomorphic there, but 
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since the numerator polynomial is of degree m{r — 1) — 1, they have simple poles over 
A = oo. To obtain the exact location of these poles and their residues, recall that the 
numerator of (1.85) was obtained by evaluating 

^■^^ (1.91) 



a(A) dP^ a{X)V, ' 

Near A = oo, 'P(A, z) is of the form: 

r J(;h^[U■^liP^-~^)+0{X)] for case (i) 

^(^' ^) = I J^[K=l(y^ - 5) + A J::=l(p^ - eu ^^^^) ^^^^^ - ^) + ^(^')] 

I for case (ii), 

(1.92) 

where the change of coordinates (1.18) has been used. Hence, for case (i), near A = oo, 

1 dz ,, Y[k^^{Pk -~Z)- X\k^x{Pk - ~Z) d~X .^ ^^ , 

-d\ ~ T- — f^f 777^ — TT —^ (1.93a) 



a{\)dP, ELiRfc^zl^fc-^) A 

and this has simple poles at oo^ (A = 0, 5 = Pi), i > 1 and ooi (A = 0, 5 = Pi) with 
residues +1 and —1, respectively. Similarly, for case (ii), 

1 dz ^^ ]lk^^{Yk-~z)-Y\^^^{Yk-~z) d\ ,^_, 

-aA ~ r — T^ 777^ — TT —, (1.93b) 



giving again simple poles at oo^ (A = 0, 5 = y^), z > 1 and ooi (A = 0, 5 = Yi) with 
residues +1 and —1. D 

It follows, since (1.82a) is essentially the Abel map, that any function on the 
Lagrangian manifold C that is symmetric in the coordinates (A^) may be expressed 
along the flow lines in terms of quotients of theta functions on the curve S. In particular, 
for the coordinates {qi{t)} themselves, we have 

Corollary 1.7. For a suitable choice of constants {ci, fi}i=2,...r, the coordinate func- 
tions {qi{t)} satsfying eq. (1.82b) are given by: 



qi{t) = In 



^(B + tU- A(oo,)-K) 



_^(B + tU- A(ooi)-K) 
where K G C^ zs the Riemann constant. 



+ e,t + /„ (1.94) 



Proof. We use the standard method underlying the reciprocity theorems relating dif- 
ferent types of abelian differentials (cf. [GH]). Namely, on the polygonization of S 
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obtained by cutting along a canonical basis {a^, b^} of cycles, we define the meromor- 
phic differential 

rfV(p) :=d{lne{A{V)-A{p)-K)), (1.95) 

where "D is the divisor Yla=lP^i formed from the g points {j>^)^=i,...g with coordinates 

(A^,^^) and p denotes the point of evaluation on S. Since d'^ has simple poles with 

residues 1 at the p^^'s, we may express the abelian sum appearing in eq. (1.82b) as an 

integral 

3 r^u r V rp ^ 

# (1.96) 

around a contour C enclosing only these singularities of the integrand, and not the ones 
at p = {ooi}, which are logarithmic branch points. Integrating by parts and deforming 
the contour to the boundary B of the polygonization of S gives 

In e{A{V) -A{p)- K)uj^ - / In 9{A{V) - A{p) - K)uj,. (1.97) 

where the Cj are small loops enclosing the poles at ooj and no other singularities. If the 
differentials iVi were normalized, the contributions to the boundary integral from the 
pairs ±a^ in B would just be constants (the discontinuity given by the theta multiplier 
over the b^ cycle), and the contributions from the b^ terms would vanish. However, 
our differentials {u)i} differ from the normalized ones by linear combinations of the 
holomorphic differentials {uJia} in eqs. (1.83). It follows from eq. (1.87) that these 
differences contribute linear terms in t with constant coefficients. The remaining terms 
in eq. (1.97) may be evaluated by taking residues at {00^)^=1^. ..r, using the results of 
Theorem 1.6 and eq. (1.87) to yield the logarithmic theta function term in eq. (1.94). 
The constants and linear terms in (1.94) are then obtained by summing those from the 
normalizing shift with those already present in eqs. (1.82b). 

Remark: The LHS of (1.82a,b) may be interpreted as an extended Abel map from S^S 
to a generalized Jacobi variety J^{S) associated to the singularized curve S obtained 
by identifying the points {00^}, where J^{<S) is a (C*)''"^ extension of J{S). The 
extended theta function for S is obtained by multiplying the ordinary theta function 
for the nonsingular curve by exponential factors in the extended directions [C]. This 
may be viewed as the source of the additional linear terms in eq. (1.94). 
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2. The Algebraic Geometry of the Symplectic Form. 

2a. The Geometric Structure of Coadjoint Orbits 

In this section, we give a geometric description of the coadjoint orbits, based on 
the results of [AHHl] (cf. also [B]). This will be done both in the reduced case, treated 
above in lb, and the unreduced case, treated in Ic. There are two parametrizations, 
corresponding to the choices Y = (case (i)) and Y = diag (Y^), Yi distinct (case 
(ii)). In the reduced case, we consider the reduced orbits Qred of Theorem 1.4, while 
in the unreduced case, we consider, for case (i), the "restricted" orbit Q^^ defined by 
eq. (1.58), and for case (ii), the full orbit Qa^o- Note that with these choices the highest 
order term of the matricial polynomial C{X) is diagonal in both cases. 

To each element of the unreduced orbit Qa/'o? we can associate a certain set of 
geometrical data. First, there is the spectral curve iS C T, defined in (1.14), which has 
the following properties (see [AHHl] for a more detailed discussion): 

Lemma 2.1. Let {zia}i=i,...n,a=i,...r-ki be the non zero spectrum of the residue ma- 
trices Ni n?=i j^ii^i ~ ^j) at X = ai. (The {Ni} are assumed to satisfy the genericity 
conditions given in Section la.) 

(1) S passes through the points (ai^Zia), is compact and lies in the linear system 
I 0{rm) \. 

(2) Generically S has an {r—ki)-fold ordinary singular point at {ai, 0); desingularising 
S at these {(ai,0)} yields a smooth curve p : S. ^' S. 

(3) In case (ii), S passes through (A, 5) = (0, YJ) (see (1.60 )). In case (i), S generi- 
cally has r distinct points over A = 0. 

A second element has already been alluded to. Let tt : T ^ Pi be the natural 
projection, where T is, as above, the total space of 0{m). We can lift the bundles 0{i) 
on Pi to T, and so to S and S. Let all these bundles be denoted 0{j) and, if F is a 
sheaf, let V ® 0{j) be denoted V{j). In a natural trivialisation over Uq C T, 0{m) 
has a basis of sections {z, 1, A, ■ ■ ■ A"^}. Given £(A), define a sheaf E, supported over 
iS, by the following exact sequence over T (cf. (1.30)) 

O^O(-m)®'^ '^^^ O®' ^E^O, (2.1) 

where /C(-2, A) = {C{X) — zl) (i.e. fx is projection to the cokernel of {C{X) — zl)). If S 
is reduced one has 
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Lemma 2.2 [AHHl]. (1) E is a torsion free sheaf over S and is generically the direct 
image of a line bundle E over S , with 

deg{E) =g + r-l (2.2a) 

(2) One has: 

H^{S,E®0{-1)) =0. (2.2b) 

In consequence, if Fx is a fibre o/tt : T ^ Pi, the restriction map 

H^{S,E)^H^{SnFx,E) (2.2c) 

is an isomorphism, and so H'^{S,E) is r- dimensional. 

Remark: For the spectral curve S of any element of Qj^^^ , we can construct the partly 
desingularised curve S mapping to S and the line bundle E over S whose direct image 
is E as follows. There is an exact sequence of bundles over P^-i: 

O^T*P^_i(l) ^^ O®'' ^ C(1)^0, (2.3) 

where (p is just the evaluation map of sections. We consider the map of sheaves over 

T xFr-i: 

T*P,_i(l)©C(-m)®" ^ C®^ (2.4) 

The cokernel of this map is supported over the union of a curve S mapping to S and 
some projective spaces over the singular points of S where the corank of IC{z, A) is 
greater than one. We take E to be the restriction of this cokernel to S. 

A third datum that can be associated to an element of QtVq is a trivialisation of 
E over A = oo. Let {ci} denote the standard basis of H^{T , O®^). If ooj is the point 
of the curve over A = oo corresponding to the j'-th eigenvector, then from (2.1), since 
the leading order term of C{X) is diagonal, /u(ei) is nonzero over oOj only when i = j. 
Then n(ei){oOi) defines a trivialisation r of the fibre of E at ooi, i.e. over A = oo. 

We can reobtain C{X) from the triple (iS, E, r) as follows. The trivialisation r, 
along with condition (2) of Lemma (2.2) allows us to fix a basis {fi} of H^{S, E) by 
the condition: 

fiiooj) = dij. (2.5) 

Then C{X) is defined as the endomorphism of H'^{S,E) (expressed in the basis {fi}) 
defined for each A by the diagram: 

H^{S,E) > H^{SnFx,E) 



£(A) 



X z 



(2.6) 



H^{S,E) > H^{SnFx,E) 
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where xz denotes multiplication by the fibre coordinate z. This is equivalent to building 
a resolution: 

Q^[0{-m)®'^ ^H^{S,E)®0{-m)\ ^^^ [O®"- ^ H'^iS.E) ® O] ^ E ^ Q. 

(2.7) 
Let U be the variety of equivalence classes of triples (5, E, r) such that 

(1) S satisfies the conditions (l)-(3) of Lemma 2.1, is generic in the sense of conditions 
(2), (3) and is such that the curve S obtained by desingularising at the (a^, 0) is 
smooth. 

(2) E satisfies conditions (1) and (2) of Lemma 2.2, and is generic in the sense of 
condition (1). 

(3) r is a trivialisation of E over the r points of S at infinity. 

Theorem 2.3 [AHHl]. There is a biholomorphic equivalence of a non empty Zariski 
open subset of Q_^f^^ in case (i) and of Q^ in case (ii) withU. 

The geometric picture then consists of a 2(^ + r — 1) dimensional space U of triples 
(5, E, r) that projects to a (^ + r — 1) dimensional space W of curves S: 

aiU^W. (2.8) 

From the discussion of the previous section, it follows that the fibres of this map form a 
Lagrangian foliation of W, since the coefficients of the spectral curve provide a (g+r—l)- 
parameter family of commuting Hamiltonians defining this foliation. The infinitesimal 
aspects of this picture are described by the following theorem. 

Theorem 2.4. (1) The tangent space at S to the space W of curves is H^{S, K^{1)), 

where Kg is the canonical bundle of S. 

(2) Fixing S, the tangent space to the fibre a~^{S) is H^{S, 0{—l)). 

Proof. (1) First order variations of a curve S immersed in T are in one to one corre- 
spondence with sections of the normal bundle A^ of 5 in T. These variations are not, 
however, entirely free: they are constrained to vanish at {A = ctiji^i^...^, and in case 
(ii) at A = oo. In other words, they must vanish at the zeroes of a section of 0{m + 1), 
and so the permissible variations of S correspond to H^{S, N ®0{—m — 1)). However, 
0{—m — 2) is isomorphic to the canonical bundle Kq- of the surface T. Using the 
notation of eq. (1.18), this can be seen from the fact that 

dzAd~X = -X-'^-^dz A dX. (2.9) 
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From the adjunction formula, we have K^ — N ® K-r, and hence the tangent space to 
W at 5 is Ts{W) = H^{S, Kg{l)). 

(2) Line bundles E on S are classified by H^{S,0*). Using the exponential exact 
sequence, 

O^Z^C^C*^0, (2.10) 

we can express infinitesimal variations of E as elements of H^{S, O). If we assume now 
a fixed trivialisation of i? at A = oo, then in considering an infinitesimal variation of 
an equivalence class represented by a cocycle (3 over Uq nUi, we can only modify (3 
by functions on Uq, and functions on Ui which vanish at A = oo. This means in effect 
that such variations are represented by the cohomology group H^{S, 0{—l)). D 

We can write down explicit formulae as follows. The image S of S lies in the linear 
system | 0{rm) \ in TPi. If S were smooth, sections of the normal bundle of S would 
be sections of 0{rm) and hence represented by polynomials f{X,z) of degree (rm), 
with the convention deg(A) = 1, deg{z) = m. The normal vector field corresponding to 
/ would be: 

/(^-^) A (2 11) 

dP/dz{X,z) dz' ^ ■ ^ 

The fact that S may be singular at (ai,0) forces us to require that / vanish to an 
appropriate degree at (Q;i,0), in order that (2.11) be finite along each branch of S 
at (ai,0), and thus represent a section of A^^. (See the discussion in Section Id.) 
Furthermore, since we want S to remain fixed at A = a^ and, in case (ii), at A = oo, 
we must require that (2.11) vanish at A = a^ and, in case (ii), at A = oo. Thus (2.11) 
is divisible by a(A). The explicit identification 0{—m — 2) ~ Kq- used here then tells 
us that the 1— form with a pole at A = oo corresponding to (2.11) is: 

^^^'^^ dX. (2.12) 



a{X)dP/dz{X,z) 



At a point (5, E, r) of W, let /3(A, z) be a cocycle on Uq fl Ui representing a class 
in H^{S,0{—1)) and let T{X,z) be a transition function for E over Uq fl Ui with 
respect to trivialisations over Uq, Ui compatible with the trivialisation r. We can write 
down a one-parameter family (Es^Ts) deforming {Eq^tq) := {E,t), with derivative at 
s = equal to the class [[3] , by choosing for Eg the transition function (with respect to 
trivialisations compatible with Ts) 

T{X,z,s)=T{X,z)e'f^^^''\ (2.13) 
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A geometric parametrization of the reduced orbits can be constructed in a similar 
fashion. In reducing, the following extra constraints must be imposed on the spectrum 
of C{X), and so on the spectral curve S: 

(3') In case (i), S intersects infinity (A = 0) at the points z = Pi , where Pi are 
distinct constants. In case (ii), the branches of S at infinity have expansions z — 
Yi + XPi + 0(A2), Pi constants. 

We then quotient by the action of the diagonal group. This action can be thought of in 
terms of changing the trivialisation r, so we define lA^eA to be the variety of equivalence 
classes of pairs (5, E), where: 

(1) The curve 5 in T satisfies conditions (1) and (2) of Lemma 2.1 and condition (3'), 
is generic in the sense of conditions (2) and is such that the curve S obtained by 
desingularising at the (ai,0) is smooth. 

(2) The torsion free sheaf E satisfies conditions (1) and (2) of Lemma 2.2 and is 
generic in the sense of condition (1) 

Theorem 2.5. There is a biholomorphic equivalence of a non-empty Zariski open 

subset of Qred with the set Wred ■ 

Again, we have a submersion crj-ed '■ ^red -^ VVj-ed onto a space W^ed of curves, with 
^red, VVj-ed of dimensions 2g, g respectively. The fibre of this map at 5 is a Lagrangian 
submanifold, and is in essence the complement of the theta-divisor in the Jacobian of 
S. 

For the infinitesimal picture, we can repeat the reasoning of Theorem 2.4. Since 
we are imposing one extra order of constraint on the curve at A = oo, the first order 
variations of the spectral curve now correspond to sections of N ^ 0{—m — 2) = Kg- 
On the other hand, we are just considering infinitesimal variations of the bundles E 
since there is no trivialization r fixed at infinity, and so we obtain 

Theorem 2.6. The tangent space at S to W^ed is H^{S,K£). Fixing S, the tangent 
space of a~J^{S) is given by H^{S, O). 

2b. Symplectic and Algebraic Geometry 

The maps U -^ W, U-^ed -^ VVred, the theorems in the preceeding sections, and 
Serre duality give us exact sequences, at a point p of W or U^ed'- 

^H^{S, Kg{l)Y -^ TpU -^ H^{S, K^{1)) -^ (2.14a) 

O^H^iS^K^y ^TpU^ed^ H^{S,Ks) ^0. (2.14b) 
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Given splittings of these sequences the tangent spaces can be written as sums of vector 
spaces of the form A* © A, where A = H^{S, K^(l)) in the unreduced case and A — 
H^{Si, K£) in the reduced case. A natural skew form u can then be defined on such a 
sum: 

uj{{a,v){h,u)) = a{u)-h{v). (2.15) 

A splitting of (2.14a,b) amounts to finding some way, infinitesimally, of fixing the line 
bundle (and trivialisation) while varying the curve. There is a natural geometric way of 
doing this. At a point (5, E, r) of U (resp. (5, E) of Wred), we extend the bundle E to 
a neighbourhood of S in T. Similarily, we extend the trivialisation to a neighbourhood 
of iSnFoo in Foo. This gives us a "background" {E, r) (resp. E) to restrict to variations 
of the curve. These extensions, and the splittings of (2.14a,b) that they define, are not 
unique. They do, however, define the same 2-forms via (2.15). (This will appear as a 
direct consequence of the proofs of Theorems 2.7 and 2.8 below). For the time being, let 
us suppose that some arbitrary choice has been made. Since these forms involve Serre 
duality, we denote them by oos and uJs,red- Let uJorb^ ^red denote the Kostant-Kirillov 
forms on the orbits as in Section 1. 

Theorem 2.7. Under the identification between U-^eA cind a Zariski open set in the 
reduced orbit Qred given in Theorem 2.5, we have 

t^red = ^S,red- (2.16) 

Proof. Let {S{x,t),E{x,t)) be a two parameter family in our space, with S{x,t) de- 
fined by the equation P(A, z, x, t) = 0, and the family of sections of E{x, t) given over 
{t/i}i=o,i by the functions {si(A, 2, x, t)}i=o,i- We shall now evaluate the symplectic 
2-form on the tangent vectors -^, -^ to U at (x,t) = (0,0). To do this, we use the 
description given by Theorem 1.4. If ^^(A^,, z,y) is the divisor of s away from A = cxd, 
so that the points (Aj^, Zi,) are given by the simultaneous vanishing of / and sq, then 
the reduced Kostant-Kirillov form evaluated on ^, -^ is 

""-' U' ^J ^ V — ^<>^) — ■ 

where the subscripts x,t denote differentiation. Now define new variables Aj^, Zi, by 

A,(0,0) = A,(0,0) 
P(A,(x,t),2(A,(x,t)),0,0) = 
SQ{Xi,{x,t),z{X^{x,t)),x,t) = 

5,(0,0) = 2,(0,0) 
P(A,(0,0),5,(x,t),a;,t) = 0. (2.18) 
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Here A;^(x,t) is the A-coordinate of the point cut out on 5(0,0) by the equation 
so(A, 2, x,t) = and Zj, represents the variation of the ^-coordinate of S{x,t) over 
A = A;y(0,0). Imphcit differentiation of the equations (2.18) and of the corresponding 
equations for Aj,, z^, aUows one to show that at (x, t) = (0, 0), 

{K)t{zi,)r, - {Xu)x{zi,)t = {Xu)t{zu)x - {K)x{zu)t- (2-19) 

Let d denote exterior differentiation along the curve 5(x, t), with (x, t) fixed. We have, 
for (x, t) smaU, 

A^ = (b A rf In So (2.20) 

27r2 Jc, 

for some suitable contour C,,, and so 



ly: 



{K)t or X = TT" f Hd In So)t or X = TT" f i^^ ^o)t or xdX. (2-21) 

27rz Jc^ 2m Jc^ 



At (x,t) = (0,0) 



^'^^'°'^" P,(A.,^.,0,0) ^ ^ 



Let 



c/A (In so)t -[x^t]. (2.23) 

«(A)Pz / J 

Then at (x, t) = (0, 0), we have 

i,X„)t{ziy)x - {\u)x{zu)t ^ J_ f p ,^ 24^ 

a(A,) 27rz Jc^ - ^ - ' 

Choose a base point Aq and cut open the Riemann surface 5(0, 0) into a 4^-gon (g — 
genus(iS)) in the standard fashion. Then, 



""^(l-l;) = i?£/ 



^ ' ' -f -E/ ^l- (2.25) 



^TTZ \ ./edge of 4g— gon • ^D 



where the Dj are contours around the points ooj over A = oo. Since the expression F 
in (2.23) is defined on the curve itself along the cut locus, the contributions of the two 
sides of the cut to (2.25) cancel, so the integral along the edge is zero. Also, the 1-form 
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^.^yp dX is holomorphic at A = oo. As above, we then write the transition function 
for the hne bundle E{x,t) to first order in (x,t) as T{X,z) e*^*+^^==, so that 

so(A, z, X, t) = si(A, z, X, t) T(A, z) 6*'^*+'^^== (1 + r(A, z, x, t)V{X, z, x, t)) (2.26) 

for some function r on Uq HUi, where sq, Si represent sections along the curve. Sub- 
stituting into (2.25) gives 

-- (!• ^) - 2^ i: [i (M) d'' ^ <'"">' ^ '"<^ ^ '■^'')) - <^ " "■ 

(2.27) 
Of the three terms in the integrand, the second gives zero since (Insi)t is holomorphic 
at A = oo, and the third vanishes after antisymmetrization {x ^^ t). Thus: 



To complete the proof, we note that: 

^) oihT ' a(x)v ^^^ simply the elements of H^{S, Kg) corresponding to the varia- 
tions of the curves in the x, t directions, and (3t, (3^^ are representative cocycles for 
the elements of H^{S, O) representing the variations of the line bundles E{x,t), 
as in Theorem 2.4. 

2) The sum of the contour integrals around the Dj is the explicit representation of 
the Serre duality pairing, when the cocycles are chosen with respect to the Uq, Ui 
covering. D 

For the unreduced case, the corresponding result is: 

Theorem 2.8. Under the identification between U and a Zariski open set in the orbit 
QNo for case (ii), and the symplectic submanifold Qj^o for case (i) given in Theorem 
2.3, we have 

^ / > I 1 V^ dPi A dPj ^n 

Case (i) c^orblQO^ =^g+9Z^ r, _T, over Q^^ 

^ ^^J ^' ^^ (2.29a) 

Case (ii) uiorb = ^s- (2.29b) 



Proof. Repeating verbatim the proof of Theorem 2.7, we have, instead of (2.27) 

EdXa A dza f d d \ 1 v^ / VxdX /^ ,, , ,, , ^nn \ / x 

(2.30) 



36 ADAMS, HARNAD AND HURTUBISE 

Again, the first term corresponds to ug, and the third term disappears after antisym- 
metrization. The second term, however, does not, since now the 1-form ?|xy can have 
a pole at A = oo. The residue of this pole at the point ooi is precisely —{Pi)x- 

Now let us recall the explicit rebuilding of C{X) from the sections of E (see 
[AHHl]). If fj is the basis of sections of E, normalised as in (2.5), then for each 
A, we can evaluate fj in the Uq trivialisation at the r points (A,2j(A)) of S above A, 
and set: 

(V'(A)).,=/,(A,z„(A)). (2.31) 

Then 

£(A)=V'-'- diag(^,(A))- V- (2.32) 

Near A = oo, '0 can be expanded: 

(V').j=5.j+A7,,+0(A2). (2.33) 

Dividing (2.32) by A"^, we have, in case (i) 

£(A) ■ A-"+i = Lo + ALi + ■ ■ ■ 

= diag(PO + A([diag(P0,7] + diag((Li),,)) +0(A2) 



(2.34a) 



and, in case (ii) 

n 

£(A)A-" = Y + A(Lo - J] a,Y) + ■ ■ • 

= diag(y,) + A([diag(y,), 7] + diag \ P^ 'Y. ^J^^ I ) + ^(^')- 

V ^=1 / (2.34b) 

The section s = (sq, si) used in our calculations is just /i, so that 

siCX, 5^(A)) = 5,1 + A7,i + 0(A2). (2.35) 

Thus, at A = (A = oo). 



In si(A,5i(A))j =0 (2.36) 

and for i 7^ 1, at A = 0, in case (i), using (1.63), 

(lnsi(A,5,(A)))^=(ln(Li),i),-(ln(P,-Pi))^ 

= {q^-l E HP^-PJ))t-{HP^-Pl)),, 

j^i, j>i (2.37a) 
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while in case (ii), using (1.64), 

(lnsi(A,5,(A))^ = ln((Lo),i), 

= {qi)u (2.37b) 

and similarly for (Insi)^;. Referring to Theorem 1.5, in case (i), the expression (2.30) 
becomes: 



LOorb 

\ at nx I 

i=2 






\ at nx I 

i=2 



ly(^MPj%—iPM^, (2.38a) 



and in case (ii): 

d d 



f2(^1MP^)^-iQ^)AP^)t^ 



Ot OX , 

^ i=2 

(a a \ "^ 

9t dxj ^^' ^2.38b) 

thus proving the theorem. D 

2c. Reductions to subalgebras 

Reductions to subalgebras of Ql{r)^ can be obtained by considering the fixed 
point sets of one or several involutions a on gl(r)"'" (cf. [AHP, AHHl, HHM]. This 
procedure can be used to obtain all the "classical" loop algebras. Coadjoint orbits in 
the reduced algebras correspond to unions of components of the fixed point sets on the 
orbits of the unreduced algebra. 

2c. 1. Involutions on Loop Algebras Induced by Involutions on g[(r, C). 

Let a : 0l(r, C) -^ ^[(r, C) be an involutive automorphism. We can define corre- 
sponding linear (resp. antilinear) involutions on gi{r)~^ by 

a(£)(A) =a(£(A)) (resp. (t(£(A)), (2.39) 
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depending on whether a is hnear or antihnear. The involutions a that occur in reduc- 
tions to the classical algebras fall into three types. 

Type (i) (antilinear involution) 

o-(x)=7^7"\ 77 = ±1 (2.40) 

(This gives, e.g., the reduction to 0l(r, ffi)"*".) The corresponding a induces an antiholo- 
morphic involution i on the "spectral surface" T: 

i(^,A) = (z,A). (2.41) 

A fixed point of a has its spectrum fixed by i. From the exact sequence (2.12) defining 
E, we have: 



(zl-£(A)) 



-^ 0{-m)®'^ ' ' " ) C'®" ^ i*E — > 0, (2.42a) 

whereas, the sheaf cf[E) corresponding to d{C) is given by: 



-^ Q^_^^®rl±}_^Wyr^^®r _^ ~(^) ^ Q_ (2.42b) 

The maps 7 : O®"^ -^ C®"", 7 : C(-m)®'' -^ 0{-m)®'' then give us, by (2.42a,b), an 
isomorphism between a{E) and i*E. If we assume that 7 preserves the diagonal form 
at infinity, we obtain (see also [HHM]): 

Theorem 2.9. For involutions (2.4O) of type (i), 

(1) The fixed point set of a on the orbit Q^„ , or on Q^ed, is a real symplectic manifold, 
with (real) symplectic form given by the restriction of the symplectic form on the 
ambient space. 

(2) Under the identification of Q_^f^^ with the set U of triplets (5, E, r) given in The- 
orem 2.3, the action of a is 



a{S,E,T) = {i{S), i*E, i*T). (2.43) 

For this case, real Darboux coordinates may also be obtained. If we choose an 
eigenvector Vq in (1.29) which is invariant under the corresponding involution on C^, 
we obtain an i-invariant divisor XK-^a^ -^a*)- "^^^ points in the divisor can be ordered 
so that 

(A2^,^2m) = (-^2^+17^2^+1)7 ^ = 1, . . ., S 

(A,,z,) = (A,,^,), u = 2s + 2,...,g (2.44) 
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for some s. It foUows that the real and imaginary parts 

1111 

(— =i?e(A2^), -j=Re{z2^), (— =/m(A2^), -J=Im{z2f,))^,=l,...s, (A^, z^)^^2s+2,...,g 

(2.45) 
are Darboux coordinates on the real submanifold of fixed points in Qred- For the 
unreduced case Q_\fQ, the remaining Darboux coordinates are similarly obtained from 

[Qi, Pi)i=2,...r 

Type (a) (linear involution) 

(t{x) = -7x^7^, 7 = ±7-^ = ±7'^ (2.46) 

(This gives, e.g., reductions to o(r, C) and sp (|,C)). Here, a induces on T the invo- 
lution: 

i(A,z) = (A,-^), (2.47) 

and so determines a corresponding map on spectral curves. For the line bundles the 
map is slightly more complicated. From the defining sequence (2.1) for E we have, for 

i*E, a{E): 

— > 0{-m)®'' (~"^~^(^^) > C)®^ — ,i*E — . (2.48a) 

-^ 0(_^)®r7(^l+^^(A))7^^^©. _^ ^^^^ _^ Q_ ^2.48b) 

Locally, we can use (2.48a,b) to represent sections of i*E, a{E) by a, 6 G O®"^, respec- 
tively. If < , > denotes the standard bilinear pairing O®^ x O®^ -^ O, consider: 

<7(-zl-£(A))a,6>, (2.49) 

where '~~' denotes, as above, the classical adjoint. It is easy to check, from (2.48a,b), 
that this projects to give a pairing oi i*E with (j{E) over i{S). (Remember that 

i{S) has equation det(— ^1 — C{\)) = 0.) Since the entries of {—zl — C{X)) lie in 
if°(iS, 0{{r — 1)to), this gives a globally defined map: 

eE®a{E) >C((r-l)m). (2.50) 

When i{S) is smooth, the adjunction formula tells us that 0{{r — l)m) is Ki(^s){'^)- 

Since in this case, the adjoint matrix {—zl — C{X)) is everywhere of rank one, we can 
show directly from the sequences (2.48) that (2.49) is surjective, and so 

a(£;)~K,(5)(2)®(r£;)*. (2.51) 
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When i{S) is not smooth then, in the generic case considered, we can identify the 
canonical bundle K^,^^ of the desingularisation i{S) of i{S) with 0{{r — l)m)[—D], 
where D is a positive divisor supported by the singularities of i{S) (see [AHHl]). On 

the other hand, {—zl — C{X)) also vanishes at the singular points of i{S) in such a way 
that the image of (2.50) in 0{{r — l)m) is also C((r — 1)?ti)[— D]. Therefore again, with 
a slight abuse of notation: 

a(£;)^K,(5)(2)®(z*(£;))*. (2.52) 

Now recall that over A = oo, £ is diagonal. We assume that 7 preserves this form. The 
trivialisation of o"(-E) is then the same as that of i*E. Summing up: 

Theorem 2.10. For involutions (2.46) of type (ii), 

(1) The fixed point set of a on the orbit Q (or Q^ed) is a complex symplectic manifold. 

(2) Under the identifications of Theorem 2.3, the action of a on the triplet {S,E,t) 
is 

a{S,E,t) = {i{S), {i:'{E)Y®K^^sp),i''T). (2.53) 

Remarks: 

i) lii{S) = S, the fixed point set of the action of a on the Jacobian of 5 is a translate 
of the Pryni variety associated to i. 

ii) It is not clear geometrically what the Darboux coordinates should be in this case. 
When S is hyperelliptic (so that there is an extra involution j) such coordinates 
can be found [AvM]. 

Type (Hi) (antilinear involution) 

a{x) = -7x'^7"\ 7 = ±7"^ = ±7^ (2.54) 

(This gives, e.g., the reduction to u(p, q)). For this case, a induces on T the involution 

i{X,z) = (X,-z), (2.55) 

and so determines a map on spectral curves. Proceeding as above, we have: 

Theorem 2.11. For involution (2.54) of type (Hi), 

1) The fixed point set of a in the orbits Q (or Qred) is a real sympletic manifold, with 
(real) symplectic form given by restriction of the symplectic form on the ambient 
space. 
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2) The action of a on the triplet (5, E, r) is 

a{S,E,r) = (i{S),K^^sp)®{¥Er,¥T). (2.56) 

2c. 2 Twisted Involutions 

Given an involution a : gl(r) -^ 0^^) ^^ can define a "twisted" involution a on 
0l(r) by 

a{C){X)=a{C{-X)). (2.57) 

The case by case study of types (i), (ii) and (iii) above can be repeated. For each of 
these, the map i induced by a on T is that induced by a, composed with (A, z) hh^ 
(—A, z). With this modification of i, Theorems 2.9 , 2.10 and 2.11 again hold verbatim. 

3. Examples. 

In the following, we examine four applications of the above analysis: computation 
of Darboux coordinates on generic coadjoint orbits of s[(2) and s[(3); finite dimen- 
sional integrable systems involving isospectral flows in 51(2)+*; the cubically nonlinear 
Schrodinger equation (NLS) and the coupled 2-component nonlinear Schrodinger sys- 
tem (CNLS). Details on how these systems arise through moment map embeddings 
from a space of rank 2 or 3 perturbations of N x N matrices may be found in [AHP] . 

3a. Darboux Coordinates for sl{2)* and sl{3)* 

As a first application of the results of Section 1, we compute Darboux coordi- 
nates on generic coadjoint orbits of the algebras sl(r), r = 2,3 by viewing these as 
Poisson subspaces of the corresponding loop algebra 5l(r)"'"*. Thus, we choose n = 1 in 
eqs. (1.2), (1.3) and, without loss of generality, cti = 0. The residue matrix A^i = Lq is 
identified as an element of s[(r)* and we consider the spectral curve 

V{X,z) = det{jC{X)-zI)=0 (3.1) 

for matrices of the form 

C{X) = XY + Lo e sl{ry , r = 2,3 (3.2) 

where Y G sl{r)* is a fixed matrix with simple spectrum. (Only case (ii), with Y ^ 
will be considered, since we want coordinates on the full coadjoint orbits.) 
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The Darboux coordinates {XfiXti}^j.=i,...g,... are determined by eq. (1.29) with 
Vb = (I7O, ...)^ G C^, and the "missing" coordinates {qi, Pi}i=2,...r corresponding to 
the spectral points over A = 00 are given by eqs. (1.61), (1.64). According to the remark 
foUowing Theorem 1.4, for each eigenvalue zofY for which Vq is in the image ofY — zI, 
one of the g + r — 1 points in the spectral divisor will appear over A = 00, requiring the 
addition of a pair (qi^Pi) of "missing" coordinates to complete the Darboux system. 
In particular, if Vq is chosen as an eigenvector of Y, as in Theorem 1.5, there will be 
r — 1 such pairs associated to points at 00 and g pairs of "finite" spectral Darboux 
coordinates. 

For r = 2, the ring of Casimir invariants is generated by tr(LQ) and the generic 
orbits are 2-dimensional. Let 

The spectral curve in this case has genus g = and Vq is an eigenvector of Y; hence, 
there are no finite spectral divisor coordinates, only the pair of "missing" Darboux 
coordinates 

q2 = In-u, P2 = a, (3.4) 

(valid for u ^ 0) corresponding to the eigenvalue 5 = — 1. Alternatively, choosing Y to 

be 

1 



^=(1 OJ' P-6) 

the curve still has g = 0, but Vq is not an eigenvector of Y, so the spectral divisor 
consists of one point at finite A. Eq. (1.29) reduces to the linear system 

X + u = 

z — a = Oj (3.6) 

providing the Darboux coordinates 

Xi^-u. Ci = ^ = --. (3.7) 

Ai u 

For r = 3, the ring of Casimir invariants is generated by trLp, trLg and the generic 
orbits are 6-dimensional. Let 









1 








-1 




y=01 0, Lo= \ u a e \ . (3.8) 
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The spectral curve in this case is generically eUiptic {g = 1) and Vq is an eigenvector 
of Y. The first pair (Ai,Ci = f^) of Darboux coordinates is obtained by solving the 
system: 

{X-z + a){X + z-b)+ef = (3.9a) 

u{X + z-b) + ev = (3.9b) 

v{X-z + a)-uf = 0. (3.9c) 

The second and third of these equations imply the first, and determine the coordinates: 

X^ = -(b-a-- + ^) (3.10a) 

2 \ u V J 

zi uva + uvb — ev'^ — fu^ 

Cl = T^ = ——r^ YTm^ 3.10b 

Ai —uva + uvb — ev^ + jw^ 

when u, v and Ai are nonzero. The remaining two points of the spectral divisor, corre- 
sponding to the eigenvalues ^ = 1, — 1 of y, lie over A = oo. To complete the system, 
we must therefore add the two pairs of "missing" coordinates: 

52 = In-u P2 = a (3.11a) 

Q3 = In-y P3 = h. (3.11b) 

Alternatively, we may pick Y so that Vq does not lie in the image oiY — zl for 

any eigenvalue z of Y; e.g. 

/O 1 0\ 

y = 1 1 . (3.12) 

Vo 1 0; 

The genus is still (7 = 1, but the number of finite pairs of spectral divisor coordinates 
is now 3. These may be obtained by solving the pair of equations 

z^ + z{v — a — b) -\- X{u — e) + ab — ef + fu — av — (3.13a) 

Xz + uz + {v- b)X + ev - bu = 0, (3.13b) 

which reduces to a cubic for z, with generically distinct roots (21,22,-23)- Setting 

^^^^f + (--^-b)z. + ab-av) + fu-fe^ = f , ^ = l,2,3 (3.14) 

u — e A,- 



*i 



gives the Darboux system. 
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Continuing similarly for higher r, the ring of Casimirs for 5l{r) is generated by 
{trLo}i=2,...,r and the generic orbits are r{r — 1) -dimensional. The genus of the generic 
spectral curve is g = ^(r — 2)(r — 1) and there are, in principle, r classes of spectral 
Darboux coordinates possible, in which the number of finite coordinate pairs is between 
g and g + r — 1. In each case their determination involves the solution of a polynomial 
equation of corresponding degree. 



3b. Finite Dimensional Systems and Isospectral Flows in sl{2) 



+* 



The moment map embedding of finite dimensional integrable systems as isospec- 
tral flows in loop algebras developed in [AHP] leads, in the case s[(2)+*, to the fol- 
lowing parametrization. In eqs. (1.1-1.3), let 

Y=(^l _l^eBi{2,C) (3.15) 

and rank(A^i) = ki = 1. Then A/o(A) G gl(2)+* may be taken of the form: 

AfQ(X) = Xyil^ = XG^{A-XI)-^F, (3.16) 

^ tti - A 
1=1 

where {Fi,Gi)i^i^,,,n are the rows of a pair F,Ge M"^^ of n x 2 complex matrices 
and A = diag(Q;i) G M"^^^ is a diagonal matrix with distinct eigenvalues (ai)i=i,...n- 
Imposing the trace-free conditions tr {G J Fi) = and using the freedom of replacing 

F,y—^d,F,, G,^-^d-^G,, rf, gC-0, (3.17) 

we may take (F, G), without loss of generality, to be of the form: 

F=-L(x,y), G=-^(y,-x), (3.18) 

where x, y G C" are viewed as column vectors. (This amounts to a symplectic reduction 
with respect to the center of 0l(2, C)+*, giving flows in sl(2, C)"""*.) The reduced orbital 
symplectic form is then just 

uj = d^^ ^ dy (3.19) 

and Af{X) has the form 



a \ A / z_^i=i \-ct, l^i=i 



.2 



A/'(A) = A c \,] + -\ ^r "T' ^r'z:T l , (3-20) 



c —a ) 2 \ y^" ^i Y^*^ ^iVi 

^i=l A — a,- /-^i=l 
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where {xijyi)i=i^,,,ri are the components of (x, y). We now also impose the reahty 
conditions 

x = x, y = y, y = F, (3.21) 

to obtain flows in sl(2,M)"'"*. The Hamiltonian systems obtained by pulhng back the 
AKS ring T{si{2)*) through the map 

Jy :R'' xR"^ — > s[(2)* 
Jy: (x,y) I — >XY + XG^{A-XI)-^F^Af{X) (3.22) 

are then Poisson commutative and, with the possible addition of certain quadratic 
constraints, coincide with those studied by Moser in [M] (cf. also [AHP, AHHl]). 
(The fibres of this map are generated by the finite group of refiections {xi.yi) i — > 
{—Xij —yi) of the coordinate axes. Since the points with {xi = 0, y^ = 0) are excluded 
from the inverse image of the orbit Qj\f by the condition rank(A^i) = 1, the resulting 
ambiguity is resolved along the fiows by continuity.) 

As an illustrative example, consider the C. Neumann system [N]. This has been 
amply studied by a variety of methods in the literature [AvM, F, Kn, M, Sch, Ral]. 
We include it here to show how our general approach reduces to the familiar results 
for this case, giving a complete separation of variables in hyperellipsoidal coordinates 
and linearization via a hyperelliptic Abel map. To obtain this system, we choose the 

matrix Y to be 

-i 



and the Hamiltonian d> to be: 



Y=[o o' ' ' (^-2^) 



0(x,y) = -tr(Ar(A)2)o = i[(x^x)(y^y) + x^Ax - (x^y)^], (3.24) 

where the subscript ( )o signifies the A^ term in the Laurent expansion around A = 
for large A. To obtain the appropriate phase space, we must also add the symplectic 
constraints: 

x'^x = 1, y^x = 0, (3.25) 

defining the cotangent bundle T*S'^~^ ~ TS'^~^ C M^"^. The Neumann oscillator 
Hamiltonian is 

if(x,y) = ^[y^y + x^Ax], (3.26) 

which coincides with (/)(x, y) on the constrained manifold. The constraints (3.25) may 
be viewed as a Marsden-Weinstein reduction under the stabilizer StahiY) C s[(2, M) 
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(cf. Section lb), in which x^x = 1 defines a level set of the moment map generating 
the fiow 

(x,y)^(x,y + tx) (3.27) 

induced by the one-parameter subgroup Stab(Y) , while y^x = defines a section over 
the quotient of the level set by this flow (i.e., of the null foliation it generates). It 
follows that the if -flow of the constrained system is obtained from the (p-G.ow of the 
free system simply by orthogonal projection of the momentum y relative to x: 

(X(t), y{t))free ^ (x(t), mUnstr. := ((x(t), y(t) - (^^l^) ^(t)) (3.28) 

from the invariant manifold deflned by x^x = 1. The equations of motion for the 
unconstrained system are 

^ = (x^x)y - (x^y)x (3.29a) 

^ = -(y^y)x - Ax + (x^y)y. (3.29b) 

These are equivalent (within a quotient by the flnite group of reflections in the coordi- 
nate axes) to the Lax equation 

^ = [B,^f], (3.30a) 

where 

S = #(Ar)+=(^_^j7 ^^jj^y (3.30b) 

The invariant spectral curve is thus given by the characteristic equation 

det (:^ - CI)) = 0, (3.31) 

which, deflning as in Section 1 

n 

^:=a(A)C, a(A):= J](A-a,), (3.32) 

determines a genus g = n — 1 hyperelliptic curve defined by (cf. eqs. (1.15), (1.22): 

z"^ -a{X)V{X) = 0, (3.33) 
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where 



= P„-iA"-i + P„-2A"-2 + ■ ■ ■ + Po (3.34a) 



1 1 

i = --5^/, = --x^x (3.34b) 

i=l 
1 " 1 

-5^a,/, = -(/> (3.34c) 



4 



and 



are the Devaney-Uhlenbeck invariants (cf. [M]). 

Applying the prescription of Section lb, with Vq = (1,0)^, the reduction with 
respect to the stabihzer of Y gives rise to the constraints (3.11) as discussed above, and 
the solutions to eq. (1.34) give us the Darboux coordinates {Xfj,XtJ.)n=i,...n-i defined 
by 

^^^ n;:;(A-A,.) ^^3^^^ 

^^ X — ai a(\) 

1=1 



2 ^ A^ - tti y a(A^) 

Thus, the spectral divisor coordinates here are just the usual hyperellipsoidal coor- 
dinates (A^), together with their conjugate momenta (C^)- The Liouville generating 
function on the isospectral foliation thus becomes 



"-1 ^-^ rXu Ip^^^ 



I I, A. I I, A. p ^ 

^ = Z^ C^ldXfJ,\p^=cst. = > , / 

11=1 u=l"'0 



a(X) 



dX, (3.37) 



and the canonically conjugate coordinates undergoing linear flow are 
dS 1"^ f^^ X^dX 

where, for our Hamiltonian (f) = 2P^_2, 

bn-2 = 2, bj =0, j <n-2. (3.39) 
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This reproduces the famihar hnearization via the hypereUiptic Abel map obtained 
through the classical methods of Jacobi (cf. [M]). 

The other classical systems treated in [M] as isospectral flows of rank 2 perturba- 
tions of a flxed matrix A, such as geodesic flow on hyperellipsoids or the Rosochatius 
system, follow identically (cf. also [GHHW, AHP, Ra2]). In all these cases, the 
spectral divisor Darboux coordinates (A^,^^) will coincide with the usual hyperellip- 
soidal coordinates, or some complexiflcation thereof, the curves will be hypereUiptic 
and the spectral invariants will be an analogue of the Devaney-Uhlenbeck invariants 
(3.21) encountered in this case. 

3c. The NLS System 



The NLS equation has two distinct forms 

+ \/—lut = 2\u\'^u 



Un 



'^xx ~r 



-lUt 



-2\uru. 



(3.40a) 
(3.40b) 



Here we shall discuss (3.40a), but (3.40b) can be dealt with in a similar fashion. The 
flnite genus quasi-periodic solutions of (3.40a) are determined by a pair of commuting 
flows in Poisson submanifolds of su(l,l)"'"* consisting of su(l, 1) -valued rational or 
polynomial functions of A. Let jC{X) be a matricial polynomial of the form 



£(A) 



a(A) 



A/'(A) 



LoA'^-i+LiA'^-2 



+ 



+ L 



n—l 7 



(3.41) 



where A/'(A) G 3l(2, C)"*"* is an element of a rational coadjoint orbit of the form (1.2), 
with y = 0, simple poles at (A = ai)i=i,...n and residues Ni of rank ki = 1. The reality 
condition implying that Af G su(l, 1)+* is equivalent to the conditions Ni, Li G su(l, 1). 
The spectral curve for this case (cf. Lemma 1.1 and Proposition 1.2) is given by 



det(£(A) - zl) =z^ + a{X)V2{X) 
where P2(A) is a polynomial of degree n — 2: 

^2 (A) = P20 + P21A + ■ ■ ■ + P2,n-2A 

Choosing Hamiltonians of the AKS type (with Y = 0): 

'a(A) 



0, 



n-2 



H^ 



H, 



A" 
a(A) 



A tr(Ar(A)2) 
A^ ix{U{\f 



-B 



2,n-3 



-B, 



2,n-4 



(3.42) 
(3.43) 

(3.44a) 
(3.44b) 
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^ives the Lax equations 

d 

ix 
d 



-£(A) = [ALo+Li,£(A)] (3.45a) 

dx 



-C{\) = [A^Lo + ALi + L2, C{\)]. (3.45b) 

dt 

The leading term Lq in C{X) is the 5u(l, 1) moment map and hence is invariant under 
all AKS flows. We choose it to be 

io = ^(j \)- (3.46a) 

Parametrizing Li, L2 as 

L,= (» I) (3.46b) 

L^ = i{Jj Is). (3.«c) 

it is easily verified that s is a constant of motion, which may be set equal to zero. On 
this level set we have 

U = u, (3.47) 

and S — I tip is constant, so we may choose 

s = (3.48a) 

S = |w|2. (3.48b) 

With these values for the matrices Lq, Li and L2, the commutativity conditions for the 
flows determined by eqs. (3.45a,b) are equivalent to the condition that u{x, t) satisfy 
eq. (3.40a), and eq. (3.45a) determines L3, . . . , L„_i in terms of w and its a;— derivatives, 
up to choices of integration constants. 

In [AHP] these flows were related to reduced canonical Hamiltonian flows via 
rank 2 isospectral perturbations of matrices as follows (cf. also [AHHl] ). Consider 
the subspace W C M"^^ x M"^^^ of the space of pairs of complex n x 2 matrices 
consisting of pairs (F, G) with columns of the form 

F = i=(z,zz), G=i=(-zz,z) (3.49) 
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where z G C" is a column vector with components (-2i)i=i,...n- The space VF is a real 
symplectic subspace of M"^^^ x M"^^^ that may be identified with C"^, with symplectic 
form 

n 

ui — —i dz^ A dz — —i NJ dzj A dzj. (3.50) 

i=i 

Fixing n real constants {ai}i=i,...n, and defining as above, A :— diag(Q;i) G M"^^", one 
constructs a moment map 

(3.51a) 



J:C"- 


--5~U(1,1)+* 




J:z- 


-^XG^{XI-A)-^F 






\ ( 7 V^ I'^-'l'' 


V" ^J' 




_ A « 1^3 = 1 \-Cj 


^3 = 1 X-aj 




~ 2 , v-n 


■ sr^n \zj\^ 




\ Z^j = l A-a, 


* Z^i=l A-Oj 



=:Ar(A). 

(3.51b) 

This is analogous to the map defined by eqs. (3.20), (3.22) above, with different reality 
conditions. The fibres are again generated by the finite group of refiections in the 
coordinate hyperplanes {z^ = 0}. Thus, via J, C"^ (minus the coordinate hyperplanes) 
provides a canonical model for an orbit Q,}^ in su(l, 1)+*. This consists of elements with 
simple poles at (ai)i=i,...n with rank ki = \ matrix residues, whose kernels and images 
are conjugate 2-vectors, null with respect to the hermitian form H{w,z) = \w\'^ — \z\'^. 
The symplectic form u) in eq. (3.50) is just the Kostant-Kirillov form cjorb on this orbit. 

The pull back under J of the AKS Hamiltonians (3.44a,b) generate commuting 
flows on the symplectic space W. The condition (3.46a) gives invariant constraints 

z^z = 1 (3.52a) 

z^z = 0, (3.52b) 

while eqs. (3.48a,b) are equivalent to 

^Az = ^aj (3.52c) 



n 
Z^ 

J = l 



Iz^ 



Az\''-z^A^z = Y,a]- 2(5^a,)2. (3.52d) 



The function u entering in eq. (3.40a) is, by eq. (3.46b) 

1 



^^--zMz. (3.53) 
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The coefficients P207 • • • 7 -^2,^-2 in tlie spectral polynomial (3.43) give n — 1 real 
functions in involution on C" ~ W. Adding the further invariant P2 — (-^0)22 = 
— I Yl^=i kjP gives a completely integrable system on W. 

The genus of the spectral curve (3.42) is g = n — 2. Since the translation term 
Y is zero, the flows on the Jacobi variety of this curve linearize the above completely 
integrable systems on C"^ reduced by the su(l, 1) action. The reduction is given in this 
case by fixing the value (3.46a) of the su(l, 1) moment map Lq, and then dividing by 
the isotropy group of Lq. The isotropy group is 5"^, whose action 



^iip 



e'^z 



(3.54) 



is generated by the Hamiltonian P2. With u given by (3.53), we see that this action 
maps u to e'^^^u and hence integration of the reduced system only determines u up to 
a phase factor. 

By considering instead the constrained submanifold, Qj^ C Q//, described for 
case (i) of Section Id, we can obtain the linearization of the flows implicitly in terms of 
complex hyperelliptic coordinates and determine u explicitly without any arbitrariness 
of phase. Since Q^ is given by setting the off-diagonal terms in Lq equal to zero, the 
two real constraints given by eq. (3.52b) give a symplectic submanifold X C C"^ (minus 
the coordinate hyperplanes) which models Q%-- The functions P20, ■ ■ ■ , P2,n-2 give a 
complete set of integrals on the 2n — 2 dimensional space X . (On this space, we have 
P2,n-2 = —P^-) The constraint (3.52a)is equivalent to choosing 

-'2,n-2 — T, 

while (3.52c,d) are then equivalent to choosing 



(3.55a) 



P 



n 



0,n-3 



i=i 



Pi 



0,n-4 






(3.55b) 



(3.55c) 



Following the prescription of Section lb, again choosing Vq = (1,0)^ and requiring 
it to be in the kernel of the matrix of co factors of {C{X) — zl), the divisor coordinates 
{A^, C^^}^^=l,...n-2 are given by 



E 

i=l 



(Xi 



a(A) 



(3.56) 
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and 

which define complex hypereUiptic coordinates {\^)^.=i,...n-2 and their canonically con- 
jugate momenta {C,ij)^=i^...n-2- 

The function u in equation (3.57) is defined in eq. (3.53) and gives the off-diagonal 
term in Li. By the results of Theorem 1.5, q2 = ln{u) must be included, along with its 
canonical conjugate P2, to complete the Darboux coordinate system on X. Thus, the 
restriction of the orbital symplectic form (3.50) to X is just 

n-2 

u = ^dX^AdC^ + dq2AdP2. (3.58) 

Notice that {X^, Ci^)u=i,...n-2, q'2, -P2 appear to give 2n— 2 complex functions on X, which 
has real dimension 2n — 2. However, since we have reduced the loop algebra to su(l, 1) 
there are reality conditions satisfied by these functions (see Section 2c, case (iii)) which 
will be preserved by the fiows. Similarly, the conditions (3.52a-d) (or equivalently, 
(3.52b), (3.55a-c)) may be imposed on the initial data, and will be preserved under the 
fiows. 

From (3.57), we see that the Liouville generating function (1.76) on the isospectral 
foliation is 



S^J2 C^.dX^\p.=cst. + 52^2 = J2j y-^h^"^^ + ^2 Inu. (3.59) 



and hence (cf. eqs. (1.79a,b)), the coordinates canonically conjugate to the conserved 
quantities ^2,0 7 • • • P2,n-2 are given by 



"-' ''^■' AVA 



1 ^ r^' x'dx 

Q2, = -J2 , n^-DfxV ^ = 0,...,n-3 (3.60a) 

^7~i-^o v-«(A)^2(A) 



n-2 „\ 



where the final term is derived using the relation P| = —P2,n-2 on X. Thus, up 
to normalization of the n — 2 holomorphic differentials appearing in eq. (3.60a), this 
represents the hypereUiptic Abel map (cf. Theorem 1.6) It follows from Hamilton's 



Q2,i = Ci, 


i < n 


Q2,n-4 = Cn-4 


-t 


Q2,n-3 = Cn-3 


— X 


Q2,n-2 = Cn-2, 
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equations for h = H^ and Ht given in (3.44a,b) that the x and t dependence of the 
QiS is given by 

1 (3.61a) 

(3.61b) 
(3.61c) 
(3.61d) 

where cq, . . ■Cn-2 are constants. Expressing the singular (3rd kind) abelian differential 
appearing in (3.60b) in terms of the Abel map and the appropriate hyperelliptic 6- 
function, we obtain, through the procedure described in Corollary 1.7, the following 
explicit formula for u{x^t) 

/ N f> /, , 6'(A(oo2,p) +tU + a;V-K) 

u{x^t) = exp(,,) ^ K eMbx + ^^) ^(A(ooi,p) + tU + xV - K) ' ^^'''^^ 

where U, V G C^ are determined as in Theorem 1.6 from the Hamiltonians h = H^;, Ht, 
respectively, 6, c are determined as in Corollary 1.7, and the integration constants 
-ftT G C, p E S are determined by initial data. This agrees with the standard ^-function 
solution of NLS (cf. [P]). 

3d. The CNLS System and Flows in su(l, 2)+* 

The CNLS system is a 2-component generalization of eqs. (3.40a,b). The real 
form we consider involves two complex functions u{Xjt), v{xjt) satisfying the coupled 
system of equations: 

mt + Wxx = 2tt(|w|^ + |f 1^) (3.63a) 

ivt + t'xx = 2v{\u\'^ + \vf). (3.63b) 

These can be obtained as the compatibility conditions for a pair of Lax equations 
of the type (1.1) with A/'(A) G 5u(l,2)+* of the form given in eq. (1.2) and Y = 0, 
corresponding, as above, to the Hamiltonians 



1 



a(A) 



2 



A' 
a(A) 



H.m = 7^ T^tr(Ar(A)^) (3.64a) 



Xn-2 



HtW = o THtr(Ar(A)2) . (3.64b) 



If, as above, we set: 



£(A) = ^Af{X) = LoA"-i + LiA"-2 + LsA"-^ + ■ ■ ■ + L^_i, (3.65) 

A 
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then Hamilton's equations again take the Lax form 

■^jC{X) =[XLo + Li, £(A)] (3.66a) 

ax 

-^^X) =[X''Lo + XLi + L2, £(A)], (3.66b) 

and the flows commute. If the foUowing invariant constraints are imposed: 

2 
3 



Lo = ^ I -1 I (3.67a) 



0-1 



Li= \ u (3.67b) 




-u P + I f P —Uj 



Uj: — I w P —vu I , (3.67c) 



12 

-uv 



the CNLS equations are obtained as the compatibihty conditions for the equations 
(3.66a,b). 

It is possible, similarly to the NLS equation discussed above, to obtain an intrinsic 
characterization of the orbit corresponding to residue matrices Ni of rank /cj = 1 as 
an open, dense subset of C^", viewed as a real symplectic space, (cf. [AHP, AHHl]). 
However, the approach developed in Section 1 allows us to treat all orbits of the type 
(1.3) on the same footing, regardless of the rank ki = 1, 2, 3. Only the explicit formulae 
(1.27) for the genus g of the spectral curve and (1.25a) for the degrees of the invariant 
polynomials will change. 

By Lemma 1.1 and Proposition 1.2, the invariant spectral curve is given by a 
polynomial equation of the general form 

det(£(A) - zl) = V{X, z) = Vr{X, z) + p(A, z) = 0, (3.68a) 

where 

Vr{X,z) := -z''+zA2{X)Vr2{X)+As{X)Vrs{X) (3.68b) 

defines a reference curve Sr, determined by the initial data for a particular solution of 
the CNLS system, and 

p{X,z) :=a{X) lza2{X)J2 ^2aA" + a3(A) f^PsaA" ) (3.68c) 

\ a=0 a=0 / 
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is of the form given by Proposition 1.2 for neighbouring curves, with p = bi-, cr = ^3, 
given by formulae (1.25a ) and {Pia, P3b)a=o,...p,b=o,...(7 are the Poisson commuting 
spectral invariants. It will be convenient to reparametrize the space of polynomials 
p{X, z) as follows. Set 

p 
a2(A) J2 P^aX'' = A.pA'^-^ + P2,p-iA"-3 + ■ ■ ■ + A.oA"-"-^ + (lower order) 

a=o (3.68d) 

as{X) Y, ^3aA" = P3,aA2-3 + P3,,_iA2-4 ^ . . . ^ p^^^^2n-3-a ^ (j^^^^ ^^^^^y 

a=o (3.68e) 

Note that the lower order terms are completely determined by the terms P2a, F'2a 
explicitly appearing. 

There can be singularities in general over the points with A = a^, which are 
assumed to be resolved as indicated in Section 1. Due to the normalizations (3.67a-c) 
the spectral curve of C{X) also has singularities at A = 00. In terms of the coordinates 
A = 1/A, the eigenvalues of A~"^£(A) have the expansion around A = 

zi{X) = y-(m2+m3)A"^ + --- 

52(A) = -^ + m2A3 + ... 

53(A) = -^ + m3p + --- (3.69) 

where for generic £, ?7i2 7^ ms. This gives the curve S a triple singularity (tacnode of 
order 3) at A = 00, and the desingularisation S is 3-sheeted over A = 00. Furthermore, 
£(A) is diagonalisable in a neighbourhood of A = 00. In terms of the algebro-geometric 
constructions of Section 2, this means that the sheaf E of (2.1) is a direct image of a 
line bundle E on S. 

The genus of the curve S is g' = p + a — 3, three less than the genus g = p + a 
given in formulae (1.26), (1.27) for the generic spectral curve in the coadjoint orbit. 
Similarily, the degree of i? is p + a — 1, three less than the generic case. Accordingly, 
we have six fewer divisor coordinates than in the generic case, and hence an insufficient 
number to provide a Darboux system on the orbit Qj\f. 

The solution to this problem is to impose six extra constraints on elements of 
the orbit, so that the dimensions of the constrained submanifold coincides with the 
number of coordinates. This must be done in an invariant way. Note that, under the 
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Lax equations (3.66a,b), the matrix £(A), A fixed, evolves by conjugation. In G/(n, C), 
the spectrum is not the only invariant under conjugation; when there are multiple 
eigenvalues, we also have the different Jordan canonical forms. For example, among 
matrices with one double eigenvalue, the generic coadjoint orbits have nondiagonal 
canonical form, and there are orbits that are two dimensions smaller, consisting of 
diagonalisable matrices. With this model in mind, we impose the following further 
invariant constraints, defining a 2(p + cr — 1) -dimensional symplectic submanifold Qs C 

(i) The spectral curve has genus three less than the generic curve in the orbit, and so 
has three extra singularities, counted with multiplicity. (3.70a) 

(a) At these extra singular points, C{\) is diagonalisable. (3.70b) 

If U denotes the space of triples (curves, sheaves, trivialisations over infinity) 
obtained from the orbit Q^ , and W is the corresponding space of curves, with the 
projection map d :U ^ W, condition (3.70a) restricts us to the inverse image under S 
of a codimension three subvariety of curves in W. For curves S in this variety, d~^{S) 
is a stratified space, with a generic stratum consisting of line bundles (+ trivialisations 
at A = oo) on the singular curve S, and other strata corresponding to direct images of 
line bundles on various desingularisations of the curve at subsets of the three singular 
points. Constraint (3.70b) then restricts us to the codimension three stratum in d~^{S) 
of direct images of line bundles on the curve S obtained from S by desingularising all 
three points. 

As remarked in section 2, it is irrelevant, for the purpose of integrating the AKS 
fiows, whether one uses the Kostant-Kirillov form cJorb or the Serre duality form ivs- Let 
(A^, Cm)m=i,--,s'' ^^ ^^^ divisor coordinates and let (g^, Pi)i=2,3 be the "extra" Darboux 
coordinates as in Section Ic. 

Theorem 3.1. The restriction to Qs of the form ujs is given, over a suitable dense 
set, by 

u;s = Yl dX^. A rfC^ + ^ % AdP, + ^Yl ^^^^ (3.71) 

/i=l i=2 i^j * -^ 

Proof. One begins by noting that Qs can be described, as in section 2, as a variety of 
generically smooth curves S, along with line bundles E (with trivialisations at A = cxd) 
defined over S. Proceeding as in Section 2, we can define a symplectic form u!s,s on this 
space, using Serre duality. As in the proof of Theorem 2.8, it follows that the explicit 
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form ofLJs,s is given by (3.71). There remains only to show that ljs,s is the restriction 

to Qs of ujs- 

This is fairly easy to see. Any line bundle E over S is the pull-back of a line bundle 
E on S. Extending this bundle to a neighbourhood of S gives a splitting of the map 
(curves, bundles, trivializations over A = oo) — > (curves), both on the constrained 
submanifold Qs and on the ambient space Q^ . Using this splitting, we now write a 
pair of vector fields {V^},=l,2 on Q%^ as {v,,e,), v, G H^{S,K{1)), e, G H\S,0{-1)), 
where {Vi} are tangent to Qs along Qs. The Vi are represented in Dolbeault cohomology 
at S by holomorphic forms, with a simple pole at infinity. Along Qg, these stay finite 
at the singularities; this is the condition of tangency. The e^, in turn are represented 
by (0, 1) forms, which vanish at infinity. On Vi, V2, "^5, ^s,s are both given by 

i^i A 62 - 'i'2 A ei, 
5 

which is obviously well behaved in a neighbourhood of Q^. D 

The constraints (3.70a,b) and the relations between P2,P^ and the coefficients of 
the polynomial p{\z) imply that the coefficients P2,p,P?,,a^ -Ps.cr-i, Pz,(j-2^ P?,,a-z 
in eq. (3.68c) can be expressed in terms of the lower coefficients {P2aiP?,b) a=o,...p-i, 
&=o,...cr-4 and P2,-P3- To apply the Liouville method, we only need to know the con- 
straints to order 1 at the CNLS curve Sr corresponding to the spectral polynomial 
Vr{\,z). The constraint (3.70a), requiring the neighboring curve to have the same 
degree of singularity as Sn, is equivalent to first order to requiring the induced section 

of the normal bundle along Sr: 

p{X,z) d 



OVr/Oz dz 
to remain finite at the singular points of Sr. 



(3.72) 



At the CNLS curve this means that, passing to the coordinates z, X, the three 
first terms in the Taylor expansion at A = of the expression 

must vanish when one substitutes z — —^. This yields the linearized constraints: 

P3,a = '-P2,p (3.74a) 

A,.-! = {hp-i (3-74b) 

Ps,a-2 - ^A,p-2. (3.74c) 
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Setting 

Pj^:-'-+Pj, J =2, 3 (3.75) 

the linear variation of the {Pj) at the CNLS curve can be computed by evaluating 
the limits of the normal vector field (3.72) along the two branches of the curve as one 
approaches the singular point. To first order: 

P, = lim . ^(A)KA) + /3(A) 

A-o (-J,(A) + Si.(A))(-2%(A) - zt(A)) 

where (j, /c) = (2,3) or (3,2), and Zj{X) is as in (3.69). This gives 

P2 + P3 = if2{0)=lP2,p (3.77) 

and, if ?n2, ms are as in (3.69), 

% ( % 

= -^[-^(A,p-3)+A,a-3], (3.78) 

where c^, di are constants. Isolating -Ps^cr-a above, and using (3.77-78) to rewrite (3.74a- 
c) we obtain 



P2,p = -I [P2 + P3) (3.79a) 

P3,a = l{P2 + P3) (3.79b) 

P3,a-i = ^A,p-i (3.79c) 

P3,a-2 = ^A,p-2 (3.79d) 

P3,a-3 = ^ Ap-3 + i{m3P2 + msA)- (3.79e) 

This expresses the terms on the left in terms of the independent complete set of integrals 
of motion P2,0j -., P2,p-i, P3,o, ■■, -P3,(t-4, -P2, P3 (at least to first order around the refer- 
ence curve, which is all we need to apply the Liouville method). Since (X^, z^)u=i,..,g', 
{qi, Pi)i=2,3 form a Darboux coordinate system "up to constants of motion", if we define 
our generating function S{\j,^ 52, ?3, -P2a, -P36, -P2, P3) as in (1.76), then the canonically 
conjugate coordinates 

«-=j^| -K-IM ■{*=lll • (3-»«> 

Or^2a ) a=0,..,p-l l^ C-fSfc ) 6=0,..,(j-4 I. C-fi ) i=2,3 
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undergo linear flow: 

Q2a = C2a — Sa,p-lX — Sa,p-2t 
Qsb = Csb 

Q, = Q, (3.81) 

where C2a, csf,, q are constants. (Up to additive constants, H^ = — P2,p-i, Ht — 
—P2,p-2-) Evaluating the derivatives (3.82), taking (3.79a-e) into account, we obtain: 

^" hjo -3z- + MX)V,,(X) ;'g2^) 

9 p\iy ^2n-3-a+b 

^^Jo -3Z^ + A2{X)Vr2{X) 

Q2 = lnu + I + imsJ (3.82c) 

Qs = l^v + I + im2 J, (3.82d) 



where /, J are defined by 



I =? o , ,A^ ,,, dX 3.83a) 



^.=1-^0 



y /-A^ \2n-6 

J = y^ 5 , ,,,^ ,,, dX. (3.83b) 

f-^Jo -^Z^ + A2{X)Vr2{X) ^ ^ 

We can check explicitly that the integrands of (3.82a,b) form a basis for the holomorphic 
differentials on the curve S which has been desingularized over A = a^ and over A = oo. 
If ooi, 002, 003 are the three points of 5" over A = 00 corresponding to ^ = y , ^, ^ 
respectively, the integrands of / + im^J , I + im2 J have only simple poles over A = 00 
with residues (—1,1,0) and (—1,0,1) respectively at (001,002,003). Proceeding as in 
Corollary 1.7, we obtain the ^-function formulae 



, .. (XT? ( ^ , ,^^(A(oo2,p)+tU + xV-K) 

u[x,t) = exp{q2) = K2 exp(e2X + d2t)- 



v{x, t) = exp((j2) = Ks exp(e3a; + ^3^) 



e{A{ooi,p)+W + xV-K) 

e{A{oos,p)+W + xV-K) 
e{A{ooi,p)+W + xV-K) 



(3.84a) 



(3.84b) 



with U, V G C^ determined from the Hamiltonians h = H^, Ht as in Theorem 1.6, 
(ci, (ii)i=2,3 as in Corollary 1.7, and the remaining integration constants determined to 
satisfy the appropriate initial conditions. 
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